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Theory of Dielectric Relaxation for the Three-Dimensional 
Polar Rotator: Lattice Models Leading to Bimodal Loss 
Curves 


John D. Hoffman and Benjamin M. Axilrod 


The characteristics of the dielectric relaxation spectrum associated with some simple 
three-dimensional lattices consisting of polar molecules of specified shape have been investi- 
gated with the object of determining whether dielectric loss curves with two distinct maxima 
(bimodal loss curves) could be predicted on a theoretical basis for such systems. For the 
types of lattice and molecular shape considered, each dipole has one stable and four meta- 
stable orientations. 

The calculations show that under certain circumstances bimodal loss curves may arise 
for pear-shaped molecules in both the body-centered orthorhombic and tetragonal systems. 
All of the models lead to a single loss peak if the barrier system is sufficiently isotropic as in 
the body-centered cubic lattice. The changes with temperature of the shape of the loss 
curves, the static dielectric constant, and the configurational entropy are discussed. The 
abrupt changes in these properties that will occur at phase transitions due to cooperative 
interaction are also considered. Qualitative predictions concerning the expected behavior 
of bimodal loss curves for three-dimensional lattices, in which each dipole has only a single 
stable orientation, are compared with experiment. 

Using an argument based on the fact that broadened loss curves in monophase molecular 
crystals become narrower with increasing temperature, it is concluded that fluctuations of 
structure are not the principal cause of the broadening of loss curves. The present theory 
leads to the correct type of temperature dependence, and this strengthens the view that the 
origin of multiple relaxation times in molecular crystals is the anisotropy of the crystalline 


| 


| 


field. 


1. Introduction 


It is well known that for many molecular crystals 
in which the polarization decay mechanism is the 
result of the reorientation of permanent electric di- 
poles the dielectric relaxation spectrum is best repre- 
sented by a set (or distribution) of relaxation times 
rather than by a single time constant. This is indi- 
eated by the fact that when the experimental values 
of the dielectric loss factor for such a material are 
plotted in the usual manner as a function of the log- 
arithm of the measuring frequency, a bell-shaped 
loss curve is usually obtained which is considerably 
broader than that predicted on the basis of a single 
dielectric relaxation time. In some instances, a 
plot with two distinct maxima (bimodal loss curve) 
is found. 

One of the main objectives of this paper is to show 
how an extremely wide set of relaxation times capable 
of leading to a bimodal loss curve can be predicted 


; theoretically for systems that, for reasons of lattice 


symmetry and molecular shape, may reasonably be 
assumed to consist of three-dimensional polar rotators 
with five orthogonal orientations. Hitherto, cal- 
culations on the three-dimensional polar rotator 
have not been closely related to details of molecular 
shape and lattice structure, and have yielded only : 
rather narrowly spaced set of relaxation times leading 
f© a somewhat broadened loss curve [1],! or a very 
general type of information not concerned with the 
details of the shape of the loss regions [2]. Bimodal 
loss curves have been predicted for one special type 
of single-axis rotator [3]. 
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' Figures in brackets indicate the literature references at the end of this paper, 
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Another objective is to deduce from the theory the 
general characteristics to be expected of bimodal 
loss curves in molecular crystals. These qualita- 
tive predictions, which are not in general sensitive 
to a special choice of models, are compared with 
experiment. 

A quantitative comparison of theory and experi- 
ment is considered to be unwarranted at the present 
time. In the first place, for reasons of simplicity of 
treatment, we are forced to deal with the body- 
centered orthorhombic (6co) and  body-centered 
tetragonal (bct) systems, whereas the best dielectric 
relaxation data now available for three-dimensional 
polar rotators are for the face-centered orthorhombic 
(feo) and face-centered tetragonal (fet) forms of the 
hydrogen halides. Nevertheless, with certain ap- 
proximations to be mentioned, the results for the 


feo and fet lattices are practically identical with those 


obtained for the bco and bet types. The second and 
more serious complication is that the question of 
the exact position of the hydrogen atoms in the 
hydrogen halide structure is unsettled, with the re- 
sult that even the most stable dipole orientation is 
uncertain. 

Several specific three-dimensional models are dis- 
cussed. All are based on pear-shaped polar molecules 
situated in a bco or bet lattice. The time-dependent 
dielectric properties of each of these models are 
illustrated by predicting the shape of the loss curves 
for certain specified arrangements of the activation 
energy barriers hindering reorientations between the 
various possible positions (sites) that each dipole 
may occupy by rotating about its lattice point. 
In addition, an equilibrium dielectric property, 
namely, the total orientational polarizability, is 








calculated for each of the models. This provides 
information concerning the behavior of the static 
dielectric constant. Consideration is also given to 
the effect of change of temperature on the shape of 
the loss regions, the static dielectric constant, and 
the configurational entropy. 

The methods used to calculate the molecular 
relaxation times and the orientational polarizability 
associated with each relaxation time have been dis- 
cussed in considerable detail for three-dimensional 
lattice models in an earlier publication [1]. There- 
fore, after describing each model and setting down 
the differential equations pertaining to it, the relaxa- 
tion times and polarizabilities are cited without 
giving the details of the calculations. The loss curves 
were calculated from these quantities with the help 
of eq (2b) of the following section. 

The results obtained in reference [2] regarding the 
general properties of the relaxation times were used 
to check the results quoted in this paper. In addi- 
tion, the general formulas valid for orthogonal six- 
site models given in the appendix of reference [2] 
were used to check the calculation of the total 
orientational polarizability. 


2. Preliminary Relationships 


If a substance has 2 single MACTOSCO PIC dielectric 
relaxation time T, its dielectric behavior in a sinusoi- 
dal electric field is described Dy the Debve (4] 
equation 

” Ae wT 
>) ’ (1) 


lines 1+ wT" 


where ¢’’(w) is the dielectric-loss factor as a function 
of angular frequency w, and Ae is the increment of 
dielectric constant associated with the relaxation 
process. A plot of ¢’’ as a function of logy w, using 
eq (1), vields a narrow and symmetrical loss curve 
with a maximum at wf=1. As has already been 
indicated, the loss curve observed for many sub- 
stances is broader than that predicted by the Debye 
equation. For present purposes, it is convenient to 
think of such data as being represented by a super- 
position of Debye-type loss curves, each of which 
possesses its own relaxation time and magnitude. 


Therefore, (1, 5] 


é’’(w)= > * AcsoTs 2a) 
G 1+0°T; 
where Ts, with B=1,2, ..., is the set of macro- 


scopic relaxation times, and Aes is the increment of 
dielectric constant associated with each mode of 
decay. 

Calculations on molecular models of the type to 
be discussed here characteristically lead to a set of 
discrete molecular relaxation times 73, and a corre- 


2 A misprint in a general formula used in the polarizability calculations occurs 
in reference [1]. The last term in eq (11), the expression for yw - F, should be cos 
t{cos € and not cos £j cos 8. The correct expression was employed in making all 
calculations. 
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sponding set of molecular orientational polariz,| 
abilities a, where B=1,2, .... As in previoy 
papers, it is assumed that the macroscopic and mal 
cular relaxation times are approximately equal, j 
37s, and that each increment of dielectric con, 
stant is nesta Romar to the corresponding orientg. 
tional polarizability so that Aes» Aas, where K jg, i 
constant. The lack of an exact proportionality j 
each case is a result of the fact that the electrie fiel; 
at a molecule (internal field) is not identical to th 
applied field for condensed phases. With the ap 
proximations mentioned above we have the workin 
equation [1,5], ’ 


\eae KS? eee 9} 
@) kK > ee Dhe 
for estimating the details of the loss curves, |] 
cases where the quantity A2sa3;~Y3\e6=€,—e, js 
small (e, and e, are the static and high frequene 
dielectric constants, respec tive ‘ly , eq (2b IS an excel! 
lent approximation. tiene. in cases wher 
the Onsager expression for the ratio of the interna 
to the applied field is valid, there are reasons for 
believing that eq (2b) holds with sufficient accuracy 
for the purpose at hand tor any value of ¢—« 
(5,6,7). 

All of the types of loss curves mentioned here an 
in section | can be reproduced by inserting the ap- 
propriate 73; and @ into eq (2b). In particular, if 
two relaxation times (or sets of relaxation times | 
are very far apart, and the corresponding polarize! 
abilities are of comparable magnitude, a bimods| 
loss curve is generated. 

In the special case where only two active relaxation time 
exist, we have been able to deduce the following useful co 
dition for the appearance of clearly apparent bimodal | 
curves: 


()(G2)>s 


max 


In this equation rt; and +r, represent the longer and short 
relaxation times, respectively, and a@max and ami, the larger 
and smaller polarizabilities. Relaxation times and _polari 
abilities conforming with this equation will yield loss curve 
that exhibit a perceptible minimum between the two los 
regions. ‘ 


3. Pear-Shaped Molecules in a Body- 
Centered Orthorhombic Lattice (0-1, 0-2) 


3.1. Basic Model 


In this section the manner in which a structur 
consisting of pear-shaped molecules situated in 
lattice may give rise to a bimodal loss curve is col 
sidered. The unit cell of the molecular model to be 
discussed is shown in figure 1 (a). Each molecul 


in the lattice is supposed to consist of a large sphe 


roidal part on which a smaller more or less hem 


spherical protuberance is attached. The protuberanee, 


is shown in black in figure 1 (a). The dipole momen! 
uw is along the line connecting the center of the larger 
spheroid and the center of the protuberance. 
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The first step is to single out a typical molecule—it 
to consider the central one shown in 
figure 1 (a)—-and select a reasonable set of equilib- 
rium orientations and relative site energies 
as the central dipole is permitted to rotate with its 
neighbors held fixed. ‘This will determine many im- 
portant features of the crystalline field. In the 
present case, it is postulated that the protuberance 
on the central dipole may fit between the four nearest 
neighbors comprising each face of the unit cell so 
that the possible dipole orientations are orthogonal. 

In the face-centered orthorhombic structure there 
are 14 rather than 6 gaps between the nearest neighbors in 
which the protuberance might conceivably fit. The six most 
which in the directions of the nezt- 
nearest neighbors, are orthogonal to one another; thus these 
sites resemble those postulated for models 0-1 and 0-2 
3.2 and 3.3 The eight unstable sites are a result of 
the protuberance fitting into triads of nearest neighbors If 
it is assumed that these eight are sufficiently unstable 
to be ignored, the site model for the fco lattice becomes iden- 
tical to that for the bco. Under these conditions, the results 
quoted here for the bco lattice with pear-shaped molecules 
will hold for the fco lattice as well 


is convenient 


sites 


ico 


pronounced sites, are 


SCC- 


tions 


sites 


The assumption is made that there is one most 
stable orientation for the central dipole, and this 
deep site is designated ‘1.”’ In the lattice, 
where a<b<e for the unit cell, the equilibrium site 
energies in the s and 2’ directions 3 and 5) 
must be the same, and similarly, the pair of sites in 
the y and y’ directions (sites 2 and 4) must have the 
same energy. However, in this lattice the 
hearest neighbors in the ae faces are somewhat more 
compactly arranged than those in the be faces, the 
pair of sites 3, 5 is less stable than the pair 2, 4. 
lence, sites 3 and 5 are referred to as “high” sites, 
and 2 and 4 are designated as “intermediate”’ 
fig. 1 (b) It has been assumed that the dipole 
cannot point in the direction owing to the re- 
pulsion of the protuberance on the central molecule 
with that on the body-centered molecule below it; 
this idea has been borrowed from the work of Powles 
[8] on_molecules of similar shape. 


beo 


(sites 


since 


sites 
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In antiparallel dipole arrays of the type shown in 
figure 1 (a) where a~b=c, the orienting torque on the 
central dipole due to the electric field exerted by its 
neighbors will generally be small. In bee lattice 
where a=b=c this torque is actually zero. In such 
situations it is safe to ignore the effect of dipolar 
interaction on the site energies, but in strongly tetra- 
gonal or orthorhombic unit cells, or certain parallel 
dipole arrays, an appreciable contribution to the 
site energies may arise from dipolar interaction. 

Many of the equilibrium properties of the site 
model shown in figure 1 (b) can be obtained without 
introducing further information. For instance, the 
total orientational polarizability can be calculated in 
terms of yu, the absolute temperature, 7, and the 
equilibrium energy differences between the sites. 
However, the relaxation times of interest here de- 
pend on the activation energy barriers arising from 
the crystalline field that control the rate at which 
dipoles turn from one site to another. The local 
free energy barrier W,,; between two adjacent sites 
i and 7 determines the elementary process transition 
probability (unimolecular rate constant) k;,; as 
ki.,=A exp[—W,/kT], where k Boltzmann’s 
constant, and A a constant.’ In a similar manner, 
ki..=A exp[—W,,/kT)]. In turn, certain sets of 
combinations of the various transition probabilities 
prove to be related to the set of relaxation times that 
characterize the system. In order to obtain the 
required molecular relaxation times, it is therefore 
necessary to determine or assume certain properties 
of the activation barrier system in such a manner as 
to fix at least the relative values of the various 
transition probabilities. 

Certain simple considerations greatly limit the 


1S 


In assuming that each transition probability is controlled by the free energy*of 
activation, we follow the ideas expressed previously by Kauzmann [Rev. Mod. 
Phys. 14, 12 (1942)] concerning the elementary dipole orientation process in a 
dielectric In the text, the term “‘energy’’ is sometimes used to replace the more 
cumbersome “‘local_free energy,’’ especially in discussing the equilibrium proper- 
ties of a sits 








. 


possible choices of the relative transition probabili- 
ties for any model: 

(a) The only elementary process is the turning of a 
dipole to an adjacent site (single-jump hypothesis, 
S-J). The S-J hypothesis has been used in all 
previous calculations [1, 2, 3, 5). 

(b) If @ and 7 are adjacent sites of equivalent 
equilibrium energy, then ky.,—k;,; as the activation 
barrier the same for each process. This also 
follows from the principle of detailed balance. 

(c) Let the probability that a dipole will turn from 
a deep site to a high site be denoted Kaeopnign, and 
the probability for the corresponding reverse process 
be Knignoaeep. It follows from the fact that the 
activation barrier for the process deep-—>high must 
be greater than that for the process high deep that 


is 


Kacep *high* C Rnign—deep- 

(d) The transition probabilities for the turning of 
a dipole from a certain site to any of a number of 
adjacent sites that possess equivalent equilibrium 
energy must be identical for reasons of symmetry 
of the molecule and of the crystal lattice. 

Even with the above, the general orthorhombic 
model with pear-shaped molecules is excessively 
complicated and difficult to deal with from a mathe- 
matical standpoint. Hence we have introduced the 
assumption that reorientations in the z-y plane are 
much more difficult to make as elementary processes 
than other allowed transitions, i. e., k,_,<k,_. or 
k,_,; this condition merely requires W,_, to be a few 
times W,_, and W,_,.*. Then for purposes of cal- 
culation we pass to the limit and set k,_, equal to 
zero. As aresult of this assumption and subsequent 
limiting procedure, processes of the type 2- 
are forbidden, but it should be noted that a dipole 
may turn from one site to another in the z-y plane 
as the result of a succession of elementary processes 


~) 
7) 


by first entering site 1. The above-mentioned 
simplification was first introduced to render the 


problem more tractable, but further considerations, 
discussed below, that this condition 
physically not unreasonable under certain circum- 


suggest is 
stances. 

Assume that the repulsion between two pro- 
tuberances virtually in contact is more severe than 
that between the body of the molecule and a pro- 
tuberance under the same conditions. First, note 
that this assumption is consistent with the forbidding 
of site 6 as a possible orientation. Second, in turn- 
ing the protuberance on the central molecule in the 
z-y plane, observe (fig. 1 that when it i 
turned half-way to the next site it comes into close 
contact at a right angle with the protuberance on a 
nearest neighbor, thus leading to a situation much 
like that used to forbid site 6. It would, therefore, 
appear to be permissible on a trial basis to forbid 
this activated state and hence elementary processes 
in the z-y plane. Note finally that the stated 
repulsion hypothesis permits the necessary elemen- 
tary processes between the sites in the z-y plane and 
site 1. It was not found possible to justify clearly 
the forbidding of reorientations in the z-y plane on 


‘ The assumption that k.,))é 
the Sco lattice (see section 3.4). 


(fa)) Is 


ky Will also lead to bimodal loss ‘curves for 
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the basis of a simple “hard sphere”’ 
introducing special shapes for the protuberance. 


The “hard sphere’? model is useful in indicating possible 
sites. However, we find this concept is not very helpful as " 
guide for determining relative site energies and _ transition 
probabilities. For reasons of simplicity, suppose the neigh. 
bors are held fixed while the central molecule is reoriented 
Then with the simplest approximation the following is found 
first, either an orientational site is possible, or there is inter- 
ference and it is not allowed; second, all allowed sites haye 
the same energy; and third, with regard to the barrier heights 
either a transition allowed and the associated barrier 
negligible, or the transition is forbidden and the energy barrier 
is infinitely high. As a second approximation in using the 
“hard sphere” or rigid, incompressible molecule picture, it 
might be assumed that a site is not forbidden if a small overlap 
occurs, but rather that the site energies relative to the most 
stable site are indicated roughly by the degree of overlap and 
similarly, the heights of the free energy barriers that must be 
surmounted are indicated by the amount of volume overlap, 
The second approximation could be modified by permitting 
the neighbors to be displaced and to rotate slightly when the 
central molecule is reoriented, a probably more realistic rep- 


is 


resentation than that previously assumed. However, this 
relaxation of perfect crystalline structure would greatly 


complicate the problem. 

An alternative to the “Shard sphere” method of estimating 
equilibrium orientations and activation barriers is to use q 
potential law of the form u() A/r6+ B/r", where different 
values of A, B, and n are assigned to the body of the molecule 
and the protuberance. Such calculations are hampered by 


a lack of information concerning the exact repulsive lay 
(i. e., the value of n) that should be used. In addition, if 
directional attractive forces such as evidently exist in the 


hydrogen bond are present, this further complicates the 
situation. Nevertheless, it is believed that if further math- 
ematical labor is to be invested in this problem, it will prob- 
ably most profitable to employ a method that involves 
attractive and repulsive potentials which vary with distance 


be 


Within the limits imposed by items (a) through 
(d) above and the absence of elementary processes 
between sites in the z-y plane, two basic types of 
barrier system are still possible for the postulated 
arrangement of sites and equilibrium energies. The 


first type, denoted 0-1, has a smaller barrier for 
jumps of the type intermediate--deep than for 
high—-deep, while the second, 0-2, has a larger 


barrier for intermediate—>deep than for high—deep. 
Specific examples of these two types are considered 
separately below. 


3.2. Model O-1 


The simplest way to study the properties of the 
beo lattice with pear-shaped molecules for the special 
case where the barrier for the process intermediate 
deep is smaller than that for the process high—>deep 
is with the help of a specific example. The specific 
barrier system chosen to illustrate the behavior of 
this model is shown in figure 1 As has been 
noted previously for the orthorhombic lattice, the 
sites in the z-y plane must, taken pairwise, differ in 
equilibrium energy. Hence the intermediate sites 
are taken to have an energy V above the deep site, 
while the high sites are set at 5V/4 above the deep 
site. Although these assignments are arbitrary, it 
is emphasized that similar results are obtained for 
models where the equilibrium energy of the high 
sites is set in the range V to roughly 2V above the 
deep site. The model is so arranged that as V is 
taken to be zero, all barriers have the value Wo. 
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The transition probabilities for the elementary 
processes are defined as follows: 
Recep »int ke, Ae Wot V)/k7 > 
Kine *deep ke Ae~Wvkt 
> (4) 
Recw *high k Ae 9V/4 HW ki 
> ’ ‘ i" VkT 
Kitgn *deep h 4 h Le . 





As a result of the aforementioned simplification 
regarding reorientations in the z-y plane, kint 
le nt. 0. It is noted further that the 
ments given in eq (4) conform with considerations (c¢) 
and (d) mentioned in section 3.1, page In 
accord with (¢) ky >k,; and ky>ks, and in accord with 
d) the same transition probability is used to describe 
the elementary occurring between the 
paired sites of equivalent energy and the deep site, 
Reus high: Note also that ky 
the opty mage barriers are identical in the model. 
Item (b) does not apply in this particular model be- 


»>high 
assign- 
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proc CSSCS 


. ]. 
e@. g., A131 KR, as 


cause there are no adjacent sites with the same 
equilibrium energy. 

It is now a simple matter to write down the differ- 
ential equations that govern the rate processes. Let 
NM, represent the number of dipoles in 
sites 1, 2, at any time ¢t. Then, bearing in 
mind the S- J hy pothesis, which forbids jumps of the 
type 3->5 or 2-4 as an elementary process, the net 
rates at w whie h fk s enter and leave each site are 





dN, dt Ak, | kN, T kN, T kiN; if kN, T kN; | 
dN, dt kN, T kN, 
dN;/dt kyNgt+ksN, > (5) 
dN, dt kN, T kN, 
dN,/dt kyNs+ksM,. 

4 
The relaxation times may readily be obtained by 


solving the above equations using the methods given 
in previous papers and the polarizabilities > may be 
obtained by using the method outlined in reference 
(1 


The results for 0-1 are as follows: 





») 2 7 
I /hry, Oy “9 ( -_ 
: 7 (1+ 29+29")\3k7 
2/4 uw 
T | ok a : ( - .); 
Ye (1+29+-29°"*) 3k7 
2 
ko 17] 24 vq) 
(l+q+4¢ 29” q) 1—39q—2¢° q ue ; 
| (| ©) 
4(1+-2y 24 vg » 
2 
kh, +8g+ 29" —yq) 
| q by 2¢" Vv~U J 2g Sze ue ) 
411+-294+- 29°") vq $k 
| Le 
af (Se) 
1+29+ 29°)? JA\3K7 J 
In the above equations g=k,/k,=e~"*? and q It is seen that for small values of g a strongly bi- 
29+ 9° —4g" . 2g Ag” modal loss curve is produced which is characterized 
The details of the loss curves for this model are by a high-frequency peak which is larger than the 


v calculating the orienta- 
tional polarizabilities and relaxation times numeri- 
ally from ha. (6) for various values of the parameter 


j=exp | /kT), inserting these in eq (2b), and plot- 
img the vad curves. The relevant data are given in 
table 1, and the loss curves obtained from these data 
sedate in figure 2. 


low-frequency one. Both loss maxima are actually 
composed of two distinct peaks which are very close 
together with the result that, considered individually, 

each of the two loss peaks is very nearly of the De ‘bye 
type. As g increases the peaks tend to come closer 


Che polarizabilitic ven in the paper are derived on the assumption that the 
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together and finally coalesce, but even at g=1 the , 
loss curve is still somewhat broader than a Debvye- : P 
type peak. As may be seen from the expression a . 
g=exp(— V/kT), the trends just noted for increasing is ; | ’ 
g hold as well for the case of rising temperature. \ : 
To the approximation used here, the static dielec- 0.01 \ 
. . . . / \ 
tric constant is given by ¢, ~e,, + Aeyorar, Where ota / \ 
Ssa3. A simple analysis shows that o&,4; (and hence a | / \ ies 
e,) will rise with increasing temperature until g Sf \ 
reaches a value of 0.17; at this point the explicitly- 
shown 1/7’ term in the expression for 44) In eq (6) 
overcomes the increase with temperature of the term 
Se ae mn ” : : yhert 
containing g inside the brackets. Thus de,/d7’ is 
positive for 0<g<0.17 and negative for 0.17<g<1. | ; é 
It should be noted that the loss curves are strongly | ~ /\ 
bimodal only in the region 0<4g<0.06, i.e., where 12 / \ Val 
de,/dT is positive, and conversely, that only a some- Pad \ - 
° 7 ° : ry | / ( 
what broadened loss curve is obtained when de,/d7 3 / \ irve: 
has a strongly negative value, as for g~™1. \ 
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hy, gh--deep ke Ae " kT’ — — 
lt will be noted that the assignments listed above 
are consistent with the considerations regarding the 
relative values of the transition probabilities men- 

. of one ~_— , »p 9 hane of the dielectric lo ‘ es for model 0-1 

tioned in section 3.1. The differential equations are Figure 2. Shay i. , 
. . -_ ? pul ! qd 
identical to those viven for 0-1, except that eq (4) Is I 
used rather than eq (4) to define the transition prob- - siti iO tnctninaan seiial ti Aa intelli tai 
abilities. The relaxation times and polarizabilities | relaxation tim been set at 10-* se¢ each diagran , 
are found to be roe | Selden fms Alien tieanaar Tina. ere ™ 
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\ Values of the polarizabilities and relaxation times 
\ omputed with eq (8) are given in table 2. The loss 
\ irves calculated for some of the g values are plotted 

n figure 3. 

[tis seen that this model leads to strongly bimodal 
iloss curves which differ from those obtained with 
model 0-1 primarily in that the largest loss peak is 
ow on the low frequency side. Noting that q 
xp(—V/KT), it is seen that the two loss peaks con- 
verge with rising temperature and coalesce into a 4 
single broadened and asymmetric peak, which then 
becomes narrower with a further increase of tempera- 
ture. The total polarizability varies with qg in a 
manner identical to that found for 0-1, with the 
esult that the loss curves are strongly bimodal only | 
n the region where the dielectric constant is low and | 
dT positive. 
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3.4. Alternative Approximation Leading to 2 














Bimodal Loss Curves in BCO Systems | -- 
At this stage, it is of interest to point out that pan Mi 
another type of approximation concerning elemen- vy spt 
tary processes in the z-y plane will lead to bimodal site's diy 
behavior for pear-shaped molecules in the dco lattice. (HIGH) ; 
As will be seen in section 4, model T-1, bimodal loss It 
curves can be obtained for pear-shaped molecules in “a un 
the dct lattice by permitting reorientations in the k 
z-y plane to be more rapid than for all other ele- | *— — SITE 5 (HIGH) Sa 
mentary processes. It can be shown for the analo- , var 
gous beo case where reorientations in the a-Y plane | } ne 
are very rapid compared with the other possible SITE 4 1 lat 
elementary processes that bimodal loss curves are poe li | clo 
generated. Detailed calculations are not given, but ugh (FORBIDDEN) dh 
this result is readily seen for the limiting case where y nae! 
the dco lattice is assumed to be nearly tetragonal so = 
that results closely similar to those quoted for model 2" - 
T-1 are obtained. Bimodal loss curves can be ob- ot! 
tained with this approximation which have peaks of | Figure 4. Orientational site diagram for pear-shaped molecules J 
equal magnitude, or alternatively, similar to those im the bet lattice. “h 
calculated for O-1. Loss curves with the type of = ; : ; 
symmetry depicted for O-2 cannot be obtained when The rate equations are rs 
it is assumed that reorientations in the z-y plane are IN 5 , : 
more rapid than all other elementary processes. yo 1-N,+k’(N.+N,+N,4+-N; ‘i 
It is clear from the foregoing that bimodal be- dt ox 
havior can arise in the dco lattice with pear-shaped : " 
molecules under widely different circumstances. dN; (k’ +2k’’)No+tkN, +k’? (N,+N;) ™ 
While it is at present impossible to state with any dt s 
certainty which of the two approximations (very 
fast or very slow reorientations In. the I-Y plane) dN; (he? +22'’\N. EN, +k’ (N,- N,) > (10) wh 
will hold for any particular real system, it does seem dt : 
reasonable to suppose that either one or the other : hin 
should be valid for some substances. It is therefore dN (ke? +2k’)N RN, +k’ (Ng +N) .e 
considered possible on theoretical grounds that some dt ing 
real systems consisting of pear-shaped molecules in : ion 
the beo lattice will exhibit bimodal behavior, par- dN (he? 4 oh’) N. LIN hk’ (N,4 N,). | 
ticularly at low temperatures, dt ; J 
4. Pear-Shaped Molecules in a Body- The relaxation times are found to be ad 
Centered Tetragonal Lattice (T-1) — ten 
ma 
In a search for lattice models that led to bimodal 7213 EB 1 Op? mal 
loss curves, pear-shaped molecules in a bet lattice wil 
were naturally considered. The molecular model is | smi 
essentially the same as that depicted in figure 1(a) Ta gprs pe I) is j 
except that a=b<c. Following the type of argu- lars 
ment given in section 3, one obtains the site energies ‘ 
and orientations shown in figure 4 In the bef ’  . pas 
lattice, all of the sites in the «-y plane have the same the +k’ J abe 
energy and are designated “high’’ sites. In order ere 
to carry out a general treatment of this model let | The corresponding polarizabilities are hol 
k deco—niah ge Y 37] Th . - 
ie je’ | 9 dinii I ig\ 3kT) con 
* hnigh--deep—A : pol: 
" f a,=0 de 
© ntgn-nten =". J ty u ? 12 
. Because site 1 is deep, Q’ kh, but the transition rm 7 77] (ser) 
probability for reorientations between the equivalent 
high sites, k’’, may in the general mathematical ..O “0 pe ) 
model have anv value. — l+-4q 3kT J mt. 
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with g=k/k’ =e-Y/*T, where V is the equilibrium en- 
ergy difference between the deep and high sites. 
Mode 4 is inactive in the dielectric relaxation 
spectrum owing to a compensatory motion of the 
dipoles. geen 

Remembering that k’>k, and that 7, is inactive, 
it is seen from eq (11) that the active relaxation 
times are always close together except for the case 
k’’>k,k’, i. e., where the transitions in the z-y plane 
are considerably more rapid than all other elemen- 
tary processes. It is worth while to consider the 
nature of the loss curves for small values of g (ordered 
lattice with low dielectric constant) and values of g 
close to unity (disordered lattice with liquidlike 
dielectric constant). When k’’>>k,k’, and g is small, 
two widely separated loss curves of almost identical 
magnitude are obtained in contrast to the unequal 
peaks obtained for models 0-1 and 0-2. On the 
other hand, if the lattice is rotationally disordered 
so that all of the sites have nearly the same 
equilibrium energy, it would appear to be unreason- 
able to assume generally that the activation barriers 
will differ greatly. Hence, in a disordered lattice of 
the type considered, the condition k’’>>k,k’ may be 
regarded at least tentatively as being physically un- 
likely with the result that bimodal behavior is not 
ordinarily to be expected in such a system. 

It may be concluded that pear-shaped molecules 
in the bet lattice will exhibit bimodal behavior 
provided the condition k’’>>k,k’ is attained. Such a 
condition is most likely to arise in an ordered phase 
where the dielectric constant is relatively low and 
de/dT positive. Thus, whereas the bco lattice vields 
bimodal loss curve with either very rapid or very 
slow reorientations in the s-y plane, the correspond- 
ing effect can be obtained for the bet lattice only when 
transitions in the z-y plane are very rapid. 


ql 


As may be seen from the expression k,;=<A exp 
W,,kT), the quantities k, k’, and k’’ will tend 
to converge with rising temperature. Thus, the 


two loss peaks comprising a bimodal loss curve will 
tend to converge as the temperature is raised. As 
may be seen from a comparison of the quantities 
@+a; and a;, the two loss peaks in a bimodal curve 
will be of very nearly the same magnitude when gq is 
small, but as g is increased, i. e., as the temperature 
is increased, the peak associated with T2.3 becomes 
larger than that associated with 75. 

The results for the bee lattice may be obtained by 
passing to the limit g=1 for the model discussed 
above. For this highly isotropic lattice it is consid- 
ered improbable that the condition k’’>>k,k’ could 
hold generally, and it is therefore anticipated that 
bimodal loss curves will not commonly arise in such 
The behavior of the static dielectric 
constant for such a system will resemble that of a 
polar liquid, i. e., ¢, will be relatively large and 


systems, 


deJdT st rongly hnegat ive 


5. Discussion 


5.1. General 


It has been demonstrated that, given certain stated 
conditions that depend on the details of the erystal- 


411810—57 


69 


line field, bimodal loss curves can arise in systems 
consisting of polar three-dimensional rotators with 
assumed equilibrium orientations and site energies 
that are consistent with certain types of lattice 
structure and molecular shape. The basic cause of 
the set of relaxation times is the anisotropy of the 
crystalline field in which the molecules are situated. 
This anisotropy gives rise to different temperature- 
dependent transition probabilities (unimolecular rate 
constants) for reorientations between adjacent sites, 
and the equations describing the net rate at which 
dipoles enter and leave each site in terms of these 
transition probabilities then lead in a natural way to 
the set of discrete relaxation times. Only a single or 
narrowly spaced set of relaxation times appears if the 
transition probabilities are all identical as is assumed 
for a highly isotropic lattice. The relaxation times 
are generally most widely separated in a lattice with 
a high degree of anisotropy. 

It appears from these studies that rather special 
forms of the crystalline field are required to produce 
bimodal behavior. Even after postulating reason- 
able equilibrium orientations and relative site ener- 
gies, and applying the rules outlined in section 3.1 
concerning the permissible range of transition 
probabilities, it was generally found that parameters 
concerned with the nature of certain activated states 
remained, It will be recalled that at this stage it was 
noted that particular choices of these parameters led 
to the widely divergent relaxation times necessary 
for bimodal behavior. We emphasize that it was 
generally difficult to really justify these particular 
choices on the basis of lattice structure, molecular 
shape, or special types of interactions. Therefore, 
despite the fact that the above-mentioned special 
choices concerning the height (or relative height) of 
the activation energy barrier between certain sites 
did not seem to be excluded for any physical reason, 
it was impossible to state with any certainty whether 
any real system with the designated crystal structure 
and molecular shape would necessarily show bimodal 
behavior. It should be remarked that the conclusion 
that highly disordered (isotropic) lattices will not in 
general exhibit bimodal behavior rests on the validity 
of the idea that in a lattice with equivalent equilib- 
rium site energies, it is improbable that the barriers 
between the various sites will differ greatly. 

It is interesting to note that the general molecular 
explanation for bimodal loss curves in polyerystalline 
specimens is consistent with the suggestion of Cole 
and coworkers [9, 10, 11, 12] that the bimodal effect 
is associated with the dielectric anisotropy of each 
single crystallite. In their view, the polarization of 
a single crystal of a substance with a bimodal loss 
curve should probably be described by two or three 
principal dielectric constants, some of which may be 
presumed to decay after the abrupt removal of the 
electric field at a different net rate. As may be seen 
especially in earlier papers [1, 3, 5], this is exactly 
the situation that prevails for single crystals in the 
present type of theory. In the process of calculat- 
ing the orientational polarizabilities, the polarization 
for each mode of decay was resolved in the x, y, and 
> directions, and it was always found that a given 








mode was much stronger along one or two of these 
axes. This holds for all of the models leading to the 
bimodal effect investigated in this paper. 


5.2. Effect of Temperature on the Dielectric Relaxa- 
tion Spectrum 


In common with almost all of the models studied 
earlier [1, 3, 5] the relaxation times converge with 
increasing temperature. The fundamental reason 
for this is that each of the set of relaxation times is 
given by the theory in terms of elementary process 
transition probabilities defined as k,,—A exp(—W 
kT); clearly the k,, must converge with increasing 
temperature, and noting that the theory usually 
shows that the active relaxation times converge as the 
k,, converge, the result quoted above follows. As a 
natural result of the the active 
dielectric relaxation times, the two loss peaks in a 
bimodal loss curve tend to merge with rising tem- 
perature. If the loss is confined to a single but some- 
what broadened peak, the peak will tend to narrow 
further with rising temperature. It will be noticed 
for the models treated that, as a result of the special 
way in which the activation barriers were defined, it 
was possible to express the relaxation times in terms 
of the parameter g=exp(—V/kT). Noting that the 
active relaxation times for these models always con- 
verge as g goes from zero to unity, and remembering 
that V is the equilibrium energy difference, the con- 
vergence of the relaxation times with increasing 
temperature in these specific cases may be regarded 
as a direct consequence of the Maxwell-Boltzmann 
distribution law. 

An analysis of the shift with frequency of the two 
loss peaks in a bimodal loss curve as a function of 
temperature for all the models studied indicates that 
the low-frequency loss peak shifts more rapidly to- 
ward higher frequencies than the high-frequency one. 


convergence of 


5.3. Effects of Cooperative Interaction and Phase 
Trarsitions 


Cooperative interaction will in general tend to ac- 
celerate the convergence of the relaxation times with 
increasing temperature, interaction 
may be thought of as causing the barrier system to 
more uniform [3, 5, 9]. For the models 
given in this paper this effect may seen a 
qualitative basis by letting Vo tend to decrease with 
temperature. This leads to a more rapid increase 
of g than if Vo were invariant with temperature, and 
the net result is the aforementioned acceleration of 
the convergence of the relaxation times 

Cooperative interaction can lead to the existence 
of sharp phase transitions, where V drops abruptly 
and g increases from a small value to one near unity 
(fig. 5 This corresponds to the rapid equalizing 
of the barrier system at the transition temperature. 
With regard to the change of Vo with temperature, 
examples close to the present case may be found in 
the literature [13]. 

It is of interest to discuss the changes to be ex- 
pected in the loss curves, static dielectric constant, 
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transition at T, (schematic 


crystal structure, and in the configurational entropy 
above and below such a sharp transition for models 
like O-1 and 0-2 

Consider first the loss curves It is clear from th 
sharp change in g at the transition that an abrupt 
change ma take place in the shape of the loss Curves, 
This change will be especially marked if the phase 
just below the transition has a bimodal loss character 
but mav be much less so if the loss curve below thy 
transition is quite narrow. The calculations also 
indicate that an abrupt decrease is to be expected 
in the “average’’ relaxation time as the temperature 
is raised through the transition, and this will cause 
the loss curve to shift toward higher frequencies 

In the simplest situation, the phase below the 
transition for models 0-1 and 0-2 will 
that above it bee, the latter corresponding to the 
hight, disordered 
points at random in the five possible sites. As 
indicated previously, the heo phase could be bimodal 
particularly at low temperatures) but it is unlikely 
that the isotropic bee phase would exhibit 
behavior. For all the models discussed in this paper 
the configurational entropy at the absolute zero of 
temperature is zero, because each dipole has but one 
possible orientation. Above the transition where 
the phase is rotationally disordered, the configura 
tional entropy the value R In 5 entropy 
units. Hence the transition beo—>bee would involve 
a# Maximum entropy of transition of R In 5 entropy 
units. The actual value will be slightly less than 
this due to the slight premonitory 
order below the transition [note the slight increase 
of g in figure 5 below the transition]. 

A more complicated set of phase transitions 1 
possible for the beo lattice: the entropy might be 
het In such a case the 


he heo. and 


such 


assumes 


cained in the steps beo hee. 


sum of the entropies of transition would still be | 


it is to be 


~R In 5 entropy units. Furthermore 
should 


expected that, if bimodal behavior 
occur, it is most likely to be found in the ordered 
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heo phase compare models 0-1, 0-2, T-1 Bimodal 
Wide caiies to the eateue « olume process. Expeth 
ental entropies may contalt fairly lar contribution due t expansion 0 


the crystal 


lattice where the protuberance | 


increase of dis- 


the 

rapl 
at tl 
at t 
tran 
tran 
sequ 
cons 
ing ¢ 
jorm 
to t] 
atiol 
loss 
tran 
syst 
neal 
they 
orth 
tow! 


5.4. 


lt 
sion 
plete 
tion: 


eur 


certs 
dict 
whi 
The 
is 
In g 

a 
at | 
the 
relat 
unlil 
syst 
relat 
coef] 

b 
ners 
with 
high 
with 


( 
will 
case 
T-1) 

T 
inely 
solic 
atur 
it is 
tho 
app 





Ntropy 


models 


om the 
abrupt 
Curves 
phase 
racter 
Ow thy 
iS also 
pected 
‘Trature 
| cause 
ences 
rw othe 
O, and 
to the 
eTrance 
S. AS 
imodal 
nlikely 

such 
paper 
ero Ol 
ut one 
where 
figura- 
itropy 
nvolve 
itropy 
; than 


of dis 


crease 


ons 8 
rht be 
se the 


ill be | 


to be 
should 
n( lered 
modal 


Experi 


' r 
ansio 


ings curves might possibly appear for the bet phase 


I > °° 
it is sufficiently ordered above the transition, but 


h behavior is not to be expected for the bee form. 

The static dielectric constant tends to increase 
below the transition even if Vis held constant. In 
the cooperative case, it will tend to do so more 
rapidly, and will rise abruptly to a liquidlike value 
at the transition owing to the sharp increase of g 
at that temperature. The above holds for the 
transition heo where V0. gz! above the 
transition, but in the case of the more complicated 
sequence beo—>bet->bec, an abrupt rise in dielectric 


suc 


hee. 


' «onstant may take place at either transition depend- 


ng on the relative stability of the deep site in the bet 
form. 

The remarks in the above paragraphs pertaining 
tothe behavior of the dielectric constant, configur- 
ational entropy, crvstal structure, and shape of the 
loss curves in the vicinity of cooperative phase 
transitions apply to the corresponding face-centered 
systems provided the eight sites due to the triads of 
nearest-neighbors are sufficiently unstable that 
they may be ignored in comparison with the five 
orthogonal sites where the central dipole may point 


3.1). 


sO 


toward next-nearest neighbors (see section 


5.4. Qualitative Predictions and Comparison With 
Experiment 


It is considered that the work with three-dimen- 
sional models with one deep site Is sufficiently com- 
plete and reneral to cite certain qualitative predic- 
concerning the nature of the bimodal loss 
eurves which can in such systems. With 
certain qualifications which will be noted, the pre- 
dictions mentioned below hold for all such models 
which have been studied (many unpublished). 
The basic reasoning leading to these predictions 
is mainly to be found in the preceding discussion. 
In general, the theory indicates that: 

a) Bimodal loss curves are most apt to be found 
at low temperatures in anisotropic phases (such as 
the orthorhombic) where the dielectric constant is 
relatively low and de, dT is positive. They are 
uilikely to arise in highly disordered and isotropic 
systems, especially where the dielectric constant is 
relatively high and has the negative temperature 
coefficient typical of a polar liquid. 

b) The two loss peaks will tend to converge with 
nereasing temperature, and the loss peak connected 
with the longer relaxation time will shift toward 
higher frequencies more rapidly than that associated 
with the shorter relaxation time. 

c) The relative magnitude of the two loss peaks 
will change with temperature, but in a few special 


tions 


arise 


eases this effect may be quite small (see model 
T-1). 

These predictions may easily be extended to 
include the changes expected in the vicinity of 


solid-state phase transitions. If the low-temper- 
ature phase is bimodal just below the transition, 
it is commonly to be expected that only a single 
though perhaps highly broadened) loss peak will 
appear in the more isotropic high-temperature 


form, and that the mean relaxation time of the 
latter will be as short or shorter than that obtained 
for either of the peaks below the transition. If the 
two peaks of the bimodal loss curve have already 
coalesced prior to the transition, a further narrowing 
may take place at the transition, but the most 
marked effect in this case will probably be an abrupt 
shift of the loss maximum to a higher frequency. It 
should be noted, however, that the persistence of 
the bimodal effect above a phase transition is not 
rigidly excluded by the theory: a possible example 
would the bet transition where the dct 
phase remained ordered. 

Some aspects of the general validity of the theory 
may be examined by comparing the predictions 
given above with what has been found experimentally 
for the feo, fet, and fee structures in the hydrogen 
halides. It is believed that such a qualitative com- 
parison is permissible, even though the calculations 
on which the predictions are based refer directly 
only to the beo, bet, and bee structures.’ First, it will 
be recalled that under certain circumstances (section 
3.1) the results for the body- and face-centered lat- 
tices may be identical. Second, the predictions listed 
are just those that are insensitive to the particular 
three-dimensional single deep site model chosen, and 
we have been unable to uncover any reasonable 
models of the type mentioned above where these 
predictions failed generally. Finally, even if the 
hydrogen halides are hydrogen-bonded, so that the 
protuberance (hydrogen) points toward a nearest 
neighbor, thus forming sheets of zigzag chains [14] 
where the sites may be nonorthogonal, there seems 
to be no valid reason for abandoning the orthogonal 
orientational model with one deep site as a useful 
first approximation. In the hydrogen-bonded “site” 
model, the deep site would simply be the primary 
hydrogen-bonded position. 

A summary of the dielectric investigations of 
Cole and coworkers [9, 10, 11, 12] on the solid 


hvdrogen halides is given in table 3. Insofar as 


be beo 


the necessary data can be obtsined from their 
work, the general qualitative predictions given 


above appear to hold. Note that bimodal behavior 
is found in the most anisotropic (orthorhombic) 
phases in accord with (a), and that the two loss 
peaks undergo changes with temperature of the 
type mentioned in (b) and (ce). 

As mentioned in section 4, bimodal behavior is 
not to be expected for a tetragonal lattice that is 
disordered to such an extent that the static dielectric 
constant behaves in a manner similar to that of a 
liquid (e. g., «, is large and de, d7T is negative). 
The tetragonel forms mentioned in table 3 exhibit 
negative values of de,/¢7, and insofar as the nature 
of the loss curves has been elucidated, only single 
loss peaks are found. It should be mentioned, 
structure similar in appearance to that depicted in figure 
1 (a) except that the dipoles are all parallel to one another [see W. J. Dulmage 
und W Lipscomb, Acta Cryst. 4, 330 (1951)]. Fortunately, the same model 
fig. 1 (b)) holds for both the parallel and antiparallel cases. At higher tempera- 
tures, HCN has a fet structure. The parallel dipole dct structure leads to the 
site model shown in figure 4, so that the theory given for model T-1 may well 
ipply to this phase. Hence, a direct application of the theory to body-centered 
systems may be possible when sufficient dielectric data are available for such 


systems. An application should not be attempted if the parallel case leads to 
ferroelectric behavior, but this is not anticipated in the case cited. 


HCN possesses a beo 
t 








TABLE 3. 


Diele ctric he havior of the hydroge nm halide 8 


and cowol ke rs 


afte r ( ‘ole 


Substance Low-temperature fco form (phase III Intermediate-temperature fct form (phase High-temperature fee form (phage ] 
II 
Bimodal at low temperatures The two loss | This phase absent in HCl; feo form con- Dielectric loss not investigated, pj. 
| peaks converge with rising temperature with verts directly to fcc form at transition electric constant similar to that of 
1 change in relative magnitude, coalescing i polar liquid (e, large, de./dT nega. 
HCI just prior to the first-order transit on I'he Live 
: low-frequency peak shifts toward higher fre 
quencies with rising temperature more rapidly 
than the high-frequency one e. low, de/dT 
positive 
Bimodal at low temperatures Relative magni- Debye-type on low-frequency side of loss Dielectric loss not investigated, pj. | 
tude of the two loss peaks change with tem curve, but high-frequency data insuffici electric constant similar to that ofa 
perature. May be bimodal up to transition ent to fully determine shape. Loss max- polar liquid («, large, de/dT nega- 
HBr , DBr Low-frequency peak shifts toward higher imum is at higher frequency than in fco tive 
frequencies more rapidly than high-frequency form. e¢, unusual wd high near A-transi 
ones. « unusually high near A-transition, but tion, de,/d T negati 
de,/dT positive 
| Loss curve highly broadened but not bimodal | Debye-type over considerable frequency Dielectric loss not investigated. Dj. 
Hi", DI in temperature range studied. «. low, de/dT range but shape of loss curve not fully electric constant similar to that of 
| positive defined on high-frequency side Loss i polar liquid («, large, de./dT nega- 
maximum is at higher frequency than tive 
in feo form Dielectric constant liquid- 
like («, large, de./d T negative 
however, that any such qualitative comparison 5.5. Bimodal Loss Curves Due to Other Causes 


may well be invalid in the immediate vicinity of 
the lambda transition between the feo and fet phases 
in HBr and DBr where the dielectric constant 
reaches very high values. These high dielectric 
constant values indicate a strong dipolar correlation 
effect not taken into account in the polarizability 
calculations. Such an effect not occur near 
the first-order transition in hydrogen chloride, and 
is small near the lambda transition in hydrogen 
iodide. The substitution of deuterium for hydrogen 
sometimes leads to large increases in the observed 
relaxation times [11;. Such a phenomenon can be 
understood qualitatively if it is assumed that 
quantum tunnelling affects some of the transition 
probabilities 

It is worth pointing out that the properties of feo 
hvdrogen chloride [9] are similar to those of model 
0-2. ‘The loss curves are of approximately the type 
depicted in figure 3, and the static dielectric constant 
rises with increasing te mperacture. Also, the entropy 
of the first-order transition has been listed by Powles 
as ?} In 4.7 entropy units, which is consistent with 
that suggested for model 0-2 

The general type of theory proposed here and in 
earlier work may be checked further by examining 
the properties of substances that possess loss curves 
with only a single maximum. The theory indicates 
that the relaxation times will generally converge with 
increasing temperature, and this leads to the predic- 
tion that broadened loss curves confined to a single 
maximum should become narrower with increasing 
temperature.’ In the cases known to the authors 


does 





(section 5.3) 


where sufficient data are available to test it. this 
prediction is generally borne out by experiment; 
there appear to be no cases where a loss curve 
broadens with rising temperature. 

* Results of dielectric investigations are frequ | ported i 
Argand diagram (Cole and Cok reular a Int ! 
tation the prediction is that if the center of the circle deserib the experimental 
point sd presse 1 below the e’ ordinate. it ould r t vard ita 
ture israised, Stated in another way, ifthe Cole and Cole parameter a 


greater than zero, it should decrease with rising temper 
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It seems necessary to mention that bimodal loss 
curves in real materials can arise from a trivial cause. 
namely, that two distinct phases, each with its own 
relaxation time, may coexist during the measure- 
ment. Under equilibrium conditions such a two- 
phase system will, of occur except at 
melting or first-order transition temperatures, but it 
is well known that two phases (at least one being 
metastable) have a tendency to coexist over a wide 
range of temperatures in certain compounds. It is 


course, not 


therefore clear that if bimodal loss curves are found | 


experimentally, it should not be assumed that they 
are a property belonging to a single phase without 
conducting proper checks aimed at detecting a second 
phase. X-ray diffraction experiments can 
distinguish between monophase and two-phase sys- 
tems. Also, if at any given temperature the relative 
magnitudes of two loss peaks depend markedly on 
the thermal mechanical history of the sample, it 
is quite likely that two more phases may be 
present, There is no indication that an of the 
experimental bimodal loss curves mentioned in this 
paper are due to the presence of two phases. 
Resonance absorption due to oscillatory motions 
of dipoles may occur at microwave higher fre- 
quencies in polar solids. If a material shows absorp- 
tion due to this effect, and further exhibits a typical 
dipole absorption at lower frequencies, a bimodal 
loss curve will result. The present theory ex hypoth- 
est deals only with dipole orientation, and does not 
account for bimodal loss curves where a high fre- 
quency peak is due to the resonance effect. Reso 
nance absorption will lead to a loss peak that 3 
narrower than that predicted by the Debye equation, 
and is restricted on theoretical grounds to frequencies 
pi ene to relaxation times which are ~10™ 
r shorter. The validity of a site model, as well 
as th a of the S-J hypothesis, depends on the presence 
of barriers which are in excess of the prev — 
thermal energy kT. To be on the safe side, 1 Is 
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best to assume that each of the Wy is at least 5kT. 
With the expression r= (h/kT)exp(W4/kKT) one then 
obtains the rough estimate that the present theory 
is valid for systems with relaxation times of ~107"° 
sec or longer. Another way of estimating the range 
the shortest relaxation time which one might expect 
due to dipole rotation must certainly be somewhat 
longer than that calculated for free rotation [5]: 
this again leads to the suggestion that the theory 
should be valid for relaxation times which are ~107"° 
sec or longer. ‘Thus, an incorrect application of the 
theorv resulting from the intrusion of resonance ab- 
sorption can generally be avoided simply by dealing 
with data obtained at frequencies not much in excess 
of 10° or 10° eyeles per second. In many cases it 
should be possible to establish the absence of reso- 
nance absorption and thus treat data obtained at 
somewhat higher frequencies. 


5.6. Fluctuations 


It has been assumed in the present type of theory 
that, as the selected central dipole undergoes re- 
orientations, the nearest neighbors remain fixed. 
This is considered to be an excellent approximation 
as long as the temperature is low and most of the 
molecules are in the ground state, and such a state 
of affairs should prevail for values of the parameter 
g=exp (—\/kT) that range from zero to roughly 
0.01 or 0.1. However, when one of the nearest 
neighbors to the central molecule is itself out of 
position owing to thermal agitation, it is no longer 
strictly correct to maintain the same set of transition 
probabilities for processes involving jumps past this 
particular neighbor. The approximation implicitly 
used here and in earlier papers in dealing with this 
situation is equivalent to assuming that such fluctua- 
tions of structure can be “averaged” for the central 
dipole, which is reminiscent of the Bragg-Williams 
approximation commonly used in equilibrium order- 
disorder [15] and orientational disorder [13] theory. 
The effect of fluctuations of this type would be most 
marked in disordered crystals at high temperatures; 
this condition corresponds to g@1 where almost all 
of the molecules are in “‘wrong’”’ orientational posi- 
ions. 

There are good reasons for believing that fluetua- 
tions of structure of the type described above lead 
toonly second-order effects in the dielectric relaxation 
spectrum of monophase molecular crystals, and that 
it is therefore misleading to attribute to this cause 
alone the broadening of loss curves in such systems. 

In the first place, it is clear that if fluctuations of 
structure played the predominant role in causing 
the broadening of loss curves, rising temperature 
would an inereased broadening of the 





cause loss 
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curves. As is brought out in section 5.4, this is 
contrary to all the presently available experimental 
evidence concerning the temperature-dependence of 
the shape of dielectric loss curves. Second, if it is 
assumed that fluctuations are the sole cause of the 
broadening of loss curves, it is very difficult to see 
how bimodal loss curves could arise in monophase 
crystals. A third point is that a few crystals ex- 
hibit a Debye-type loss curve (single relaxation time) 
at all temperatures. Such a phenomenon can be 
understood, at least to a first approximation, in 
terms of an isotropic crystalline field (or certain 
special site models), but adopting the fluctuation 
viewpoint, it would appear to be necessary to arbi- 
trarily postulate the total absence of fluctuations at 
all temperatures to explain such a result. 

Noting that the present theory leads in a natural 
way to the correct type of temperature-dependence 
of the shape of the loss curves, and that it can yield 
both bimodal and Debye-type curves, it is our 
conclusion that it is generally more satisfactory to 
attribute the existence of multiple dielectric relaxa- 
tion times in monophase molecular crystals to the 
anisotropy of the crystalline field in which the 
dipoles reorient. 


The authors thank F. Buckley, 5. G. Weissberg, 
J. 1. Lauritzen, and A. D. Franklin for making a 
number of helpful suggestions during the course of 
this work. 
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Classification of Perovskite and Other ABO.-Type 
Compounds’ 
Robert S. Roth 


\ classification of A*2B*4O; compounds has been made on the basis of ionic radii of the 


constituent lons, 


A graph of this type for the perovskite compounds can be divided into 


orthorhombic, pseudocubic, and cubic fields, with an area of ferroelectric and antiferroelectric 


compounds superimposed on the 


cubic field. 


The structures of solid solutions between 


various perovskite compounds cannot be completely correlated on the basis of this simple 


two-dimensional chart. 
pounds has been attempted, using 
plotted as the third dimension. 


to compounds of the type A*?B**Q;, which crystallize with perovskite-type structures. 


Therefore, a new type of classification of the perovskite-type com- 
a three-dimensional graph with polarizability of ions 
This three-dimensional graph has been applied successfully 


In 


addition, a classification of double oxides of trivalent ions has shown thateall compounds of 
this type having the perovskite structure can be divided into two groups, with rhombohedral 


and orthorhombic symmetry. 
been found. 


1. Introduction 


A partial survey of the reactions occurring in 
binary oxide mixtures of the types AO: BO, and 
4,0,:B,0, has been conducted as part of a program 
of fundamental research on ceramic dielectrics. In 
addition, binary, ternary, and quaternary reactions 
between end members of the AO-BO, compounds 
have been studied. Combinations of oxides 
giving the simple formula type ABO, were selected 
because of current interest in ferroelectric ceramics. 
Previous publications have covered the data com- 
piled on many of these systems, such as CaO-TiO, 
1|.° PbTiO.-PbZrO (2], BeO with Zr( Vo, TiOg, and 
CeO, [3], MgO-CaO-SnO, and TiO, [4], MgO-ZrO,- 
TiO, and CaQ-ZrO,-TiO, [5], alkaline earth oxides 
with GeO, [6], PbTiIO,-PbZrO,-PbO:SnO, — and 
PbTiO,-PbHfO, [7], alkaline earth oxides with UO, 
, and divalent ions with ZrO, [9]. This paper is 
designed primarily to present new data and a new 
method of representation of these structure types 
and to coordinate the information already published. 


also 





s 


2. Sample Preparations and Test Methods 


In general the materials used were the highest 
available purity of the component oxides, varying 
from about 98.5- to 99.9-percent purity. The start- 
ing materials, in sufficient quantities to give either 
a 10.0-¢ sample or a 1.0-¢ sample, depending on the 
availability of the raw materials, were weighed to 
the nearest milligram, mixed together, and formed 


into l-in. or 4-in.-diameter disks at a pressure of 
5,000 Ib/in.2 The pressed disks were fired for 4 hr 
at 1,100° C on platinum foil in an air atmosphere, 


ising an electrically heated furnace wound with 
S0-percent-Pt 20-percent-Rh wire. 

These disks were then ground and remixed, and 
new disks for study of the solid-state reactions, about 


:in. high, were formed at 15,000 lb/in.? in either a 


in. or 4-in.-di + mold. The specimens wert 
rin. or 4-1n.-dlameter mold. 1e specimens were 
a 

‘ This work has been sponsored s part of a program for Improvement of 
Flezoelectric Ceramics by the Office of Ordnance Research of the Department 


{the Armv 


? Figures in brackets indicate the literature references at the end of this paper. 
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No ideal cubie perovskites of the A**B*8O; compounds have 


then fired in a conventional platinum-wound quench 
furnace, or refired in the original calcining furnace. 
The quenching technique was used, whenever appli- 
cable, because it has been observed that sharper 
X-ray patterns are often obtained by very fast cool- 
ing of the specimen. It is known that phase transi- 
tions in the perovskite structures are usually com- 
pletely reversible and cannot be frozen in by quench- 
ing. Thus the perovskite structures discussed are 
the room-temperature stable forms, although other 
materials treated in this manner would generally 
contain the high-temperature forms, if any. The 
final firing temperature ranged from 1,250° to 1,550° 
C and was maintained constant for a given length of 
time. Equilibrium conditions were usually reached 
in less than 3 hr. Equilibrium was believed to have 
been reached when X-ray patterns of the specimen 
showed only a single phase or when the pattern of a 
multiple-phase specimen showed no change with suc- 
cessive heat treatment. The fired disks were exam- 
ined by X-ray diffraction, using a Geiger-counter 
diffractometer employing nickel-filtered copper.radi- 
ation. The X-ray data reported can be considered 
to be accurate to+0.001 A when three decimal 
places are recorded, and to £0.01 A° when only two 
decimal places are recorded. The results reported 
in this work are generally the data obtained at room 
temperature from specimens treated in a manner to 
form a matured ceramic dielectric body. 

In the case of specimens containing an ion that was 
apt to oxidize, a neutral atmosphere was maintained 
in the furnace. These ions were V**, U**, and Ce*; 
the neutral atmosphere used was either helium or 
argon. The furnaces used were similar to those 
previously mentioned, but modified to maintain a 
neutral atmosphere. 


3. The Perovskite Structure Type 


The structures of the perovskite-type compounds 
have been studied by many workers. The exact 
nature of the structures involved is still in doubt 
in many cases. Conflicting reports in the literature 
make difficult the job of assigning the correct struc- 








. 


ture, or even symmetry, to any particular perovskite- 
type compound. The simplest case of the perov- 
skite structure is that of a simple cubic cell with 
one ABX; formula unit per unit cell. In this case 
the A ions are at the corners of the unit cell with 
the B ion at the center and the negative ions occupy- 
ing the face-centered positons. The space group 
is O}-Pm3m and has the notation G—5 of the Struk- 
turbericht [10]. As will be seen in the following 
discussion of actual perovskite-type compounds, 
very few binary oxides have this simple cubic 
structure at room temperature, although many of 
them assume this structure at elevated temperatures. 
Many modifications of this simple structure have 
been proposed to account for the X-ray patterns 
observed for different perovskite-type compounds. 
Among these modifications a simple doubling of the 
unit cell has been used [11]. A monoclinic modi- 
fication suggested by Naray-Szabo [12] for the type 
compound perovskite (CaTiO;) has been shown to 
be actually orthorhombic by Megaw [13]. The 
latest modification of this structure was made by 
Bailey [14] and quoted by Megaw [15], using an 
orthorhombic modification with the 6 parameter 
approximately twice the pseudocubic cell edge and 
with the a and ¢ parameters approximately ,2 times 
the pseudocubic cell. It will be seen that the X-ray 
patterns of a great many compounds can be indexed 
on the basis of this type of distortion of the perov- 
skite structure. The compounds CaTiO,;, CaZrQs, 
and CaSnO, have been indexed previously on this 
basis and the d-spacings and (hkl) values published 
4, 5|. However, it should be noted that a dis- 
crepancy exists between the reported space group 
for CaTiO;, Pemn [15], and the published indexing 
of these patterns. Although CaSnO, contains no 
reflections that are not allowed by this space group, 
both CaTiO, itself and also CaZrO, show diffraction 
peaks, which can only be indexed on the basis of 
forbidden reflections. These peaks are especially 
prominent for CaZrO,; and indicate that perhaps a 
different space group is required for the CaTiO,- 
tvpe orthorhombic modifications of the perovskite 
structure. 

In addition to those types of distortions that 
require a multiplication of the pseudocubic cell, 
there are several other distortions in perovskite- 
type structures that require only a slight modifi- 
cation of one or two parameters, resulting in tetrag- 
onal, orthorhombic, and rhombohedral symmetries. 
It is these latter modifications that are of most 
interest in the study of ferroelectric forms of the 
perovskite structure, especially in the two groups 
of mixed oxides under discussion, A**?B*'O, and 
A**B*.O,. Other ferroelectric perovskite com- 
pounds, like Nat'Nb**O,, have been suggested to 
have multiple-type unit cells [16] even larger than in 
the CaTiO, type. 

A classification of the perovskite-type structures 
on the basis of the radii of the constituent metallic 
ions has been attempted by several workers {17, 
18, 16, 19]. It has been pointed out by Keith and 
Roy [19] that it is not advisable to attempt to 


classify more than one valence group of structure 
types on any one diagram. In other words, separate 
diagrams are needed for the compounds of the type 
At**BYHO,, At*B*0,, At'B**O,, etc. It seems quite 
probable that perovskite-type structures with F- 
or some other such ion, in place of O-? would also 
require separate diagrams for a reasonable classifj- 
cation. These diagrams still do not lend themselves 
to absolute classifications, as some compounds simply 
do not seem to fit. It has been suggested by Roberts 
(20, 21] that the polarization of the constituent ions 
may play an important role in determining the sym- 
metry formed by a given pervoskite structure. 4 
classic example of polarization determining crystal 
type over and above ionic radii may be found jn 
the phenacite-olivine structure types where the Zn* 
ion is found to be out of place in a radii classification 
|22]. It is this notion of polarization that is used 
in the present paper in an attempt to make a better 
classification and diagrammatic representation of the 
perovskite structures. 


4. Results and Discussion 
4.1. Mixed Oxides of Divalent and Tetravalent Ions 


a. General 

For the purpose of convenience most of the data 
accumulated has been put in the form of tables, 
although some of the more controversial composi- 
tions are also discussed in the following sections 
It is known that CdTiO, crystallizes in both the 
perovskite- and ilmenite-type structures [23]. As the 
perovskite structure requires a relatively large A ion 
and the ilmenite structure a relatively small A ion, 
it can be seen that all ABO; compounds with Ti 
and an A** ion larger than Cd*? should have the 
perovskite structure, whereas those with A*? ions 
smaller than Cd*? should have the ilmenite structure. 
In all such cases where compound formation takes 
place, this is found to be true. The values used in 
this study, for the radii of the metallic ions under 
consideration, are listed in table 1. The composition, 
structure type, and symmetry, where known, for the 
mixed oxides of divalent and tetravalent ions are 
listed in table 2. The parameters of the unit cell 
are given for all those materials studied. Informa- 
tion obtained in the present study is shown in bold- 
faced type, and the author’s comments on previous 
work are given in italics. In general, references to 
other works are made only when the _ particular 
compound or mixture has not been studied in the 
present work or when results of other investigators 
disagree with the results found in the present investi- 
gation. 

b. CeO, :TiO; 

A specimen of CeO,:TiO, heated in air has vielded 
an X-ray pattern of CeO, plus a_ perovskite-type 
phase. It seems highly unlikely that either CeQ, or 
TiO, would reduce, in air, to the point of forming an 
ABO, structure. It seems more likely that an 
A*®B*°O, structure is the phase found in the X-ray 
pattern of classification of mixed 
oxides of trivalent ions). 
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TABI B i. Radii and polarizability, where 


used.of metallic ions pertinent to this study 














Ion Radius « Polar Ion Radius * 
zat ty 
l) 
i 1 A 
Pa? 1. 34 70. 6 M1 0. 80 
rt 1. 20 91.8 F¢ 74 
Eu’ 1. If Zn 74 
Sr* 1.12 49. 1 Co 72 
Ca 1.99 1.9 Ni* 49 
Cd i; 6.9 Me 66 
Be $5 
t 
\ 0. 51 Bi 0. 96 
Ga 6 Gd 07 
Cr 63 Sm 1. 00 
I 4 Nd?* 1.04 
( st) Cet 1.07 
Ir 82 Lat 1.14 
Q2 
Petr 
c+ 0. 1¢ Snt4 0.71 
sj+4 y. Hf 78 
(; ‘ 53 Zr*4 7u 
Mn 60 Cer4 4 
Ve 63 U+ 17 
| tis I 1. 02 
s Unless otherwise indicated, radii of the ions ar« 
taken from Ahret sl 
The radius of Eut? has been adapted from the 
value given by Green [31] and corrected for radius of 
18) L40A 
Che radius of In* and S« e both given as 0.81 
4 by Ahrens [30 As ti cannot be shown on the 
liagram, In** is used as 0.82 A and Sc#8as0.89 A. This 
s thought to be justified as the compounds of Sct are 


said to have slightly smal parameters than the 
responding compounds of In*3 [19] 
Unless otherwise indicated, polarizability of the 
s taken from Roberts [20 
* The polarizability of ¢ was estimated as ex- 
plained in text 
ABLE 2. 
I itment 
| 
Is er 1 lime 
~ { nr 


bal O 
PoCO 
“co 
CaCO 


Manco 
FeCO 
ZnCO 


VMeCo 


c. SrZrO; 


The compound SrZrQ,; has generally been regarded 
as a cubic perovskite [24, 16, 25]; although Naray- 
Szabo [26] has referred it, like CaTiO;, to the mono- 
clinic system. The X-ray paitern of this compound 
shows diffraction lines other than those found in 
simple cubic structures, definitely indicating some 
type of distortion of the unit cell. In addition, the 
“cubic”? peaks in the SrZrO,; pattern have a ten- 
dency to split into two or more peaks (not Ka, and 
Ka). The Ka, doublets cannot be clearly dis- 
tinguished in the SrZrO; X-ray pattern, although 
they are quite clearly shown in patterns of the true 
cubic perovskite specimens like BaZrQ,; or SrTiO. 
The split in the pseudocubic peaks resembles those 
found in PbZrO, for a pseudotetragonal distortion 
with c/a 1, although there are a few possible dis- 
crepancies in intensity values. On the other hand, 
the distortion may be of the type found in CaTiQs. 
As observed in at least six different specimens, the 
SrZrO, ‘‘pseudocubic”’ peak at a value of h?+-4?+P= 
6 is split into a close doublet with the low-angle side 
of the doublet lower in intensity than the high-angle 


side. The same peak in PbZrO, has the intensities 
of the doublet reversed. Although CaZrO; and 


CaTiO, have four peaks very close together in this 
region, the two on the low-angle side are of lower 
intensity than the two on the high-angle side. Such 
discrepancies, noted above, make it impossible to 
assign SrZrO, to either the PbZrO, or CaTiO, struc- 
ture. Calculations of (hkl) values for both types of 
cells, as well as theoretical intensity values for each 
structure, would have to be made before a final de- 
cision could be made. 


Mixed oxides of the ly pe A 2B 4() 
Structure type Symmetry References and discussion 
Mixed oxides containing CO 
Aragonite Orthorhom bi Palache, Berman, and Frondel [39] 
do do Do 
do do Do 
do do Do 
Calcite Rhombohed: Do 
do do Do 
do do Do 
do do Do 
do do Do 
do do Do 
do do Do 
Mixed oxides containing SiO 
Unknown Present work. Specimen supplied by E 
Levin of NBS staff. 
Unknown Winchell [40] 
Pseudowollastonite Unknown Present work. Specimen supplied by E. 7 
Carlson of NBS staff 
Pseudo wollastonite Unknown Present work. Specimen supplied by H. S 


Unknown 
Enstatite 
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Orthorhombic 


Yoder of Geophyscial Laboratory, Carnegie 
Insat., Washington, D. C. 

Silverman, Morey. and Rossini [41] 

Present work. X-ray data obtained from 
specimen of Bishopsville meteorite supplied 
by H. S. Yoder. a=18.222, b=8.810, ¢ 
189 A 





TABLE 2 Vired oxides of the type A* ‘BAO Continued 





Toler Hleat treatment 
ince 
Composition factor for Structure type Symmetry Reference ind discussior 
perov lemperature met 
sk ite ( her 
structur 
Mixed oxides containing GeO 
: {1,235 l Pseudo wollastonite Uknown Roth [6 Symmetry almost heragonal 
BaGeO 1.00 - - . 
1700 l Pseudo wollastonite Unknown Roth | Heated at 20,000 psi 
PbGeO 0. 95 775 { Unknown Present work. Probably nonequilibrium 
SrGeO 2 1,235 l Pseudo wollastonite Unknown Roth | 
CaGeO XN 1,235 l Unknown Roth |6 Probably nonequilibrium 
CdGeO 87 1,100 I Unknown. Present work. Probably nonequilibri 
ZnO0:GeO 7s 1,100 0. 66 No 1:1 compound Present work. Speciment contains ZnoGeGQ 
wille mite type, plus GeOos, quartz type 
MeGeO 7 1,23 l Enstatite Orthorhombic Roth [6 a 8.661, b=8.954, « ie A 
BeO:GeO 4 No information [-} compound between BeO and GeO 
seems likel robably only BeoGeO, forms 
Mixed oxides containing MnO 
CaMno 0 84 Perov kite Cubic W ird, Gushee, MeCarrol, and Ridge 42 
| Perovskite , Cubic Yakel |43 Symmetry reported as cubie wit! 
doubled cel 
Mixed oxides containing VO 
BaVO 0. 95 1,200 (Helium 0.5 Unknown Present work. 
SrVO SS 1,200 (Helium Unknown Present work. 
CaVO 83 1,300 (Helium Perovskite Orthorhombic Present work. « 7.54 { 
. Perovskite Cubic Ward, Gushee, MeCarrol, and Ridgely [42 
MeO:VO>. 72 Unknown King and Suber [44 
BeO:VO él Unknow! King and Subs 4 
Mixed oxides containing TiO 
BaTiO 0. 93 1350 2 Perovskite Tetragonal ca >| Present work. « 189, c= 4.029 A 
PbTiO SS 1,100 l Perovskite letragonal ca >| Jaffe, Roth, and Marzullo {7 x9 
EuTio 87 Perovsk ite Cubic Brous, Fankuchen, and Banks [45 
SrTiO Sf 1,300 l Perovskite Cubic Present work. a 04 A 
CaTioO SI 1,7 l Perovskite Orthorhombic Coughanour, Roth, Marzullo. and Sennett 
i a 381, b=7.644, ¢ b4 i 
CdTio | 81 {1,265 ] Perovskite Orthorhombic Present work. «a a” Mi, 19 A 
CdTio i ' Iimenite RKhombohedra Posniak and Barth [23 
MaTio 1,250) l Ilmenite Rhombohedral Present work. «a ‘ l.a f 
FeTiO 73 Iimenite Rhombohedr Wells [24 
ZnO: Tio 73 1,24 l No 1:1 compound Present work. Specimen cont s Zm TiO, 
spinel type, plus TiQOns, rutile 
CoTioO 72 1,400 2 Ilmenite Rhombohedral Present work. \-Fay pattern too diffuse 
measure parameters 
NITiO 71 1,400 2 Iimenite Rhombohedral Present work. « {48 A, a 
MeTioO 70 1,5 2 Iimenite Rhombohedral Present work. « iA 
BeO: TiO Ho 1,800 No reaction Lang, Roth, and Fillmore |3 
\lixed o les containing SnO 
BasnO 0. 92 1.450 I Perovskite Cubic Present work. «a=/ {A 
SM) l No 1:1 compound Jaffe, Roth, and Marzullo |7 S pecime 
tains PbheSnO,, knou structure ] x 
snd) 
1650 24 No reaction Starting material PbO and SnO 
PbO:SnO go 650 ( Unknown Distorted fluorite-type structure, probably 
; PbSnO,. rting material precipitated hy- 
drated lead sta ite 
S00) PbsnO Unknow! Coffen [4 Start naterial precipitated hy- 
drated lead sta te 
Perovskite etragona Naray-Szabo [2 
Perovskite Cubic Present work. Commer specimen fired 
SrsnO & } $034 A 
. Perovskite Cubi Megaw [13], Coffeen [46 
Perovskits CaTiOs typ Naray-Szabo [26 
CaSnO St) 1,550 l Perovskite Orthorhombic Coughanour, Roth, Marzullo, and Sennett 
j a=5.518, b=7.884 64 A 
{1,375 I No compound Present work. Specime contained 0 
SnO»; CdO rolatilized. A compound of 
} CdsnO0 prot l é sts hele this tempera 
Cd8n0O a ture 
Unknown Coffeen [46 Decomposition of CdSnOQ; w 
Perovskite CaTiOg type Naray-Szabo [26 
MnoO:Sn0 . 0.74 1,250 l No reaction Present work. \in0O added as \inC Os. Spe 
men contained SnQ>o, MnyO 
ZnSnO; (? 42 Unknown Coffeen [46]. Published pattern pro! thly spinel 
plus SnO 
4 {1 300 l No 1:1 compound _.. Present work. Specimen contained CoSn0y 
CoO0:Sn0O re spinel plus SnO 
{ Coffeen [46]. Published pattern that of a mit 





ture of spinel md SnO 


78 





TABLE 2. 


Mixed oxides of the ty pe A 


2B+H40 


loler lleat treatment 
ince 
Compo on factor for Structure type Symmetry 
perov- lemperature lime 
skits ( hi 
tructure 
Mixed oxides containing SnO2— Continued 
NisnO 70 Unknown 
, SOK it Ilmenite Rhombohedral 
ALL is No 1:1 compound 
} 
(7e0 MgsSnO Hao 
{ 11,55 l No 1:1 compound 
Ged 
orms QI Unknown 
BeOQ:SnO2 oY ° 
42 Mixed oxides containing HfO,¢ 
wi 
PbH{O 0. 84 Perovskit« Pseudotetragonal* 
SrHfOs-_... 82 Perovskite CaTiOs type 
CallfO 78 Perovskite CaZrO; type 
Mixed oxides containing ZrO»! 
i 
42 . , nae 
BaZrQOs-. ). 8S 1,450 l Perosvkite - . Cubic. 
\ l Perovskite Pseudotetragonal 
PbZrO os Perovskite Pseudotetragonal 
Perovskit« Pseudotetragonal 
( l Perosvkite Pseudocubic 
} M 2 Perovskite Probably orthorhom- 
SrZrO 5 bic. 
Perovskit Cubic. 
nnett : 
Perovskite CaTiO type 
} . - 
{ i CaZrO 77 ‘) 6 Perovskite Orthorhombic 
CdO0:Zro f Perovskite CaTiOs type (ortho 
rhombic 
‘) 
MnO:ZrO 14 ; No 1:1 compound 
7nO:ZrO ( " No reaction 
CoO:ZrO Hs No reaction 
NiO:Zro ( No reaction 
ri0:ZroO 120 (Helium ; No 1:1 compound 
MgO:ZrO ‘ l 1s No 1:1 compound 
BeO: ZrO f SUM 0 No reaction 
Mixed oxides containing CeO» 
vably , . 
dh i l Perovskite Pseudocubic 
| 2 Perovskite Pseudocubic 
ih BaCeO 0.8 
Perovskite _. CaTiOg type 
Perovskit« Cubic 
fired 22 ( No reaction 
PhO: CeO Perovskite CaTiOs type 
srCeO 2 4 Perovskite Orthorhombic 
inett Perovskite CaTiOs type 
2 ] No reaction 
0 Cad -¢ No reaction 
iy i0 oO 2 Perovskite CaTiOg type 
“Ta 
No reaction 
Ys is CdO0: CeO 2 Perovskite CaTiOs type 
1 l No reaction 
pe No reaction 
MgeO:CeO 62 Perovskite CaTiO 
vel 
nO, BeO: CeO 1 1,800 0.5 No compound 
nit 
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References and discussion 


Coffeen [46]. 

‘oughanour, Roth, Marzullo, and Sennett |4). 
X-ray pattern too diffuse to measure param- 
eters 

‘oughanour, Roth, Marzullo, and Sennett [4]. 
Specimen contains MgoSnQ,, spinel, plus 
snOe 

‘oughanour, Roth, Marzullo, and Sennett {4}. 
Specimen contains Mg28nQ,, spinel, plus 
SnOs. 

Coffeen [46] 

No studies reported. Reaction unlikely from 
analogy to BeO systems. 


~ 


a 


Shirane and Pepinsky [47] 

Naray-Szabo [26]. Unlikely that pure H{O: 
was available at this time. 

Curtis, Doney, and Johnson [48]. 


Present work. a=4.192 A. 

Jaffe, Roth, and Marzullo {7}. a=4.156, 

C=4.190 A 

Sawaguchi, Maniwa, and Hoshino [49]. 

Sawaguchi, Shirane, and Takagi |50).« 

Present work. a=4.099 A. Starting material 
commercial SrZrOs, also contains ercess 
Monoclinic ZrO, plus unknown phase. 

Present work. a=45.792, b=8.189, c=5.818 A, 
Starting material SrCO; plus Monoclinic 
ZrOs, no ercess phases. 

Wells [24], Wood [16], Shirane and Hoshino 
25) 

Naray-Szabo [26]. 

Coughanour, Roth, Marzullo, and Sennett 
5). a=5.587, b 

Naray-Szabo [26]. As many of the compounds 
listed in this reference do not actually erist, 
CdZrQ; cannot be accepted as a true com- 
pound without further erperimental evidence 

Roth [9]. Contains cubic ZrO: solid solution 
plus MneQy, 

Present work. 

Dietzel and Tober [51]. 

Dietzel and Tober [51] 


SUUS, C=5.798 


Roth {9}. Contains cubic ZrO: solid solution 
plus TiO Specimen supplied by F. Brown, 


Jet Propulsion Lab., Calif. Inst. of Tech. 
Coughanour, Roth, Marzullo, and Sennett | 5) 

Contains cubic ZrO: solid solution plus MgO. 
Lang, Roth, and Fillmore {3} 


Present work. a=4.587 A 

Present work. a@=4.587 A Possibly ortho- 
rhombic.' 

Naray-Szabo (26) 

Hoffman [52], Wells [24], Wood [16]. 

Present work. 

Naray-Szabo [26]. 
not erist 

Present work. a=5.986, b=8.531, c=6.125 A. 

Naray-Szabo [26]. 

Present work. 

Keith and Roy [19] 

Naray-Szabo [26] 
not exist. 

Keith and Roy [19] 

Naray-Szabo | 26). 
not exist. 

Present work. 

Keith and Roy [19]. 

Naray-Szabo [26]. Compound probably does 
not erist. Such a compound could only be an 
antiperovskite structure. 

Lang, Roth, and Fillmore (3). 
solid solution of BeO in CeQOs>. 


Compound probably does 


Compound probably does 


Compound probably does 


Possibly a 








TABLE 2. Vixred oxides of 


the type At? BHO Continued 


Toler- Heat treatment 
ance 
Composition factor for Structure type Symmetry References and discussior 
perov- Temperature rime 
skite Cc hr 
structure* 
Mixed oxides containing UO 
Balto 0.82 | 1,900 (Argon 0.5 Perovskite Pseudocubic Lang, Knudsen, Fillmore, and Roth [8 
a { 387 A Similar to BaCeO vith rery 
small splitting of diffraction peaks 
srUO 75 1,900 (Argon Perovskite Orthorhombic Lang, Knudsen, Fillmore, and Roth {8} 
a=6.01, b=8.60, « f 7 
{1,800 (Argon Perovskite Orthorhombic Lang, Knudsen, Fillmore, and Roth [8 
Cato 71 a=5.78, bh=8.29, « 5.97 
Rare earth Cubie Alberman, Blakely, and Anderson [27] 
Mg0O:U Op» 61 1,800 (Argon 5 No reaction Lang, Knudsen, Fillmore, and Roth |s 
Bed: UO 52 1,800 (Argon No reaction Lang, Knudsen, Fillmore, and Roth |8) 
Mixed oxides containing ThO»:* 
BaThoO 0. 80 : Perovskit Cubic Hoffman [52], Wells [24], Wood [16 
Perovsk ite CaTiOs typ Naray-Szabo [2t 
» As the radii given by Ahrens [30] are closer to Pauling radii than to Goldschmidt radii, the values listed for the tolerance factor are, in general, slightly smaller 
than those calculated by other authors (Keith and Roy [19], Wood [16]), using modified Goldschmidt radii 
Probably only CaO can form a perovskite compound with MnO2 (Ward, Gushee, McCarrol, and Ridgely [42 It seems unlikely that CaMn0Os has a simple 
cubic structure, and it is tentatively left on the border between orthorhombic and pseudocubic types 
rhis structure is related to cubic perovskite as follows: a + 2a’, )™2a’, c™+/2a 
1 Due to the lanthanide contraction, Hf** has a radius only slightly smaller than Zr*4, and the compounds containing Hf{O2 are probably mostly isostructural wit 
the corresponding ZrO»: compounds 
lhe pseudotetragonal structure with c/a<1 is actually orthorhombic with a= +2a’, b=2 2a’, and ¢=2c’ 
No compounds of the ilmenite structure type have been found with ZrO, or, indeed, with any B** ion larger than Sn** 
« Stronge anomaly of the dielectric constai:t at about 230° C, antiferroelectric below this temperature 
rhe exact nature of the distortion has not been determined, as the amount of distortion is small here are definite splits in the diffraction lines as well as excess 
peaks other than those for a simple cubic structure 
No work has been reported on UO: systems with NiO, CoO, Zn O, MnO, CdO, or FeO 
Published X-ray pattern is that of a mixture of UO» solid solution and perovskite compound Pattern published for Cael O, is actually the perovskite-typ 
compound 
* Perovskites of the following oxides with ThO: have been listed by Naray-Szabo [26]: BaO, PbO, SrO, CaO, CdO, and MgO. Of these reported compounds 
probably only BaThOs really exists (Keith and Roy [19 
Although no ThO: compounds were examined in the present work, it seems likely that BaThO:; would have the same type of structure as is found in BaCeO 


ind BaUO 
d. CaUO; 


The compound CaUO, also has a CaTiO, structure. 
However, a 50:50 molecular mixture of CaO and 
UO, does not vield an equilibrium perovskite com- 
pound. The formation of the perovskite compound 
in the CaO-UO, system takes place only in the 
presence of excess CaO. This phenomena has been 
discussed by Alberman, et al. [27] and by Lang, 
et al. [8]. The perovskite compounds of the type 
A**B*‘O, have been assigned a minimum tolerance 
factor (t) of 0.77 by Keith and Roy [19]. The toler- 
ance factor for the perovskite structure, as de- 


scribed by Goldschmidt [17], is 
— Rs | R, 
y 2 ( RP» T Po) 


where 
t tolerance factor for perovskite structure, 
R,=ionie radius of larger cation, 
R»,=ionic radius of smaller cation, 
y=ionic radius of oxygen 1.40 A 


This tolerance factor for CaUQ, is equal to 0.71 and 
is the only known A**B**O, perovskite-tvpe com- 
pound with ¢ considerably less than 0.77. Some 
Ca*? may actually substitute for U** in the B posi- 
“tion, accompanied by a partial oxidation of U** to 
U*®, to give a structural formula Ca(U?7*US*,_, 


Cas?,-,))O;. Such a structural substitution results 


in a total increase in the tolerance factor. However, 
a replacement of 50 percent of U** by U*® and Ca” 
increases the tolerance factor to only 0.726. Thus 
it can be seen that the lower limit of 0.77 for the 
tolerance factor in A*t*7B**O, perovskite structures 1s 
not absolutely correct. 


e. BaGeO 


The ABO, silicate compounds cannot form perovy- 
skite-type structures, as the silicon ion is much too 
small, always occurring in tetrahedral coordination. 
The B ion must occupy an octahedral position with 
respect to oxygen, even in the most distorted struc- 
tures, in order to be classified as a perovskite. As 
the germanium ion is larger than the silicon ion and 
is known to occur in octahedral coordination, a study 
was made of the reaction of GeO, and various divalent 
metallic oxides. The only known octa- 
hedrally coordinated Ge** is one of the polymorphic 
forms of GeOQ,. GeO, can be readily formed in the 
rutile-type structure by heating the quartz-type 
polymorph at 700° C and 20,000 psi. An attempt 
was made to form a titanate-type BaGeO, (with the 
germanium ion in octahedral coordination) by heat 
ing the pseudowollastonite form to 700° € at 26,000 
psi, but no change in the crystal structure was 
observed. Differential thermal analy ses of BaGeO 


case of 


indicated a phase transformation at about 340° C. 


— 
A specimen of rut le-type GeO was prepared in jt manner by A, Van 

Valkenburg and E. Bunting of the Bureau staff ae 
Che differential thermal analysis was performed by E. Newmar Bure 


staff. 
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This phase transformation is apparently very dis- 
ruptive, as a pressed pellet of BaO:GeO, results in 
finely divided white powder of crystalline BaGeO, 
when heated to 1,000° C and furnace cooled. <A 
|-percent addition of Fe,Q, to BaGeQO, results in an 
intact ceramic specimen with an entirely different 
X-ray pattern, which may represent a stabilized 
form of high-temperature BaGeQs. 


42. Classification of ABO, Compounds Containing 
Double Oxides of Divalent and Tetravalent Ions 


classification of structure types based on the 

constituent ionic radii can be made for all the com- 
pounds of the A**B*'O, type. However, the prop- 
erties of the lid solutions formed by a mixture of 
two or more compounds cannot always be correlated 
by a simple radi classification. In an attempt to 
expls ain some of these discrepancies a third property, 
polarizability of the constituent ions, has been used. 
It is exceedingly difficult, in the present state of 
knowledge, to arrive at correct values for the polari- 
zability of a given ion in a given structure tvpe. The 
values used here have been taken mainly from the 
results published by Roberts [20, 21] and are listed 
in te ody However, as Roberts did not give a value 
for Cd**, this ion was assigned a value that would 
fit the observations of other workers. Both Pauling 
[28] and Kettelaar (29] indicate that the polariza- 
bility of the Cd** ion in a given structure type is 
between those of Bat? and of Sr**. Therefore, the 
values given by Kettelaar [29] for Bat’, Cd**, and 
Sr*? were compared, proportionately, to those given 
by Roberts for the Bat? and Sr*? ions in the peroy- 
skite structure, and a value of approximately 0.56 
A’ was assigned to the Cd** ton. 

A basic division of structure types can be made on 
a two-dimensional chart of ABO; compounds, where 
the radius of the A ion is plotted as the ordinate 
and the radius of the B** ion as the abscissa. In 
figure 1 and table 1, all values of the radii, except 
for Eu**, are taken from Ahrens [80] and represent 
the ionic size for sixfold coordination. As pointed 
out by Wood [16], a correction for coordination 
would tend to move boundary lines but not to distort 
the picture of the basic classification. The value for 
Eu** has been taken from Green [31] and corrected 
fora radius of O-? of 1.40 A. In general, a division 
made between structural types containing 
large and small <A The boundary of this 
division usually falls near a value of 1.0 A, and often 
the Ca 0.99 A) or Cd** (0.97 A) compounds occur 
in both structure tvpes. For the carbonates, all 
compounds containing divalent ions larger than Ca 
have the aragonite structure, whereas those with 
ions smaller than Cat? have the calcite structure. 
In the case of the silicates and germanates, the 
pseudowollastonite structure type occurs for most of 
the compounds containing large divalent ions, 
Whereas the enstatite structure is found for the small 
ions like Mgt? (0.66 A 

Except for the Ba and Sr vanadates described in 
table 2, all other compounds described in this study 


can be 


Ons. 
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according to the constituent ionic radit, 


, Compounds studied in the present work that have the structure shown by 
the areas bounded by dashed lines; &, compounds not studied in the present 
work that are assumed to have the structure shown by the areas bounded by 
lashed lines; @, compounds with complex or unknown structure type; X, posi- 
tion of compositions studied in the present work that do not form A#?B+4O; 
compounds; @, the presence of a compound of Fe'TiQs, and the absence of a com- 
pound at the composition ZnO:TiO: are indicated together, as Zn*? and Fet? 
have approximately the same radius, 


can be divided into perovskite and ilmenite structure 
tvpes. The boundary between the two types occurs 
at a value near that of the Cd*? ion and is shown as 
horizontal for the purpose of convenience. As no 
ilmenite compounds are known with a B** ion larger 
than the radius of Sn** (0.71 A), the lower right 
corner of figure 1 is shown as an area of no compound 
formation. In all cases the boundary lines between 
structure types are only an approximation, as not 
enough information is available on mixtures of com- 
pounds of two different structure types. In general, 
there is only limited solid solution, if any, between 
compounds of different structure types, although 
there is considerable or complete solid solution be- 
tween compounds of a given type. An exception is 
found in some mixtures of compounds on adjacent 
sides of the boundary between orthorhombic and 
cubic (or tetragonal) perovskite compounds. 

A classification of those compounds 
perovskite-type structures is shown in figure 
The At?B*O; perovskites can be divided into two 
general types, those that are cubic or have distortions 
requiring an elongation of one or more parameters, 
and those that have distortions that require some 
multiplication of the basic perovskite unit cell. 
The latter type is mainly represented by the CaTiO, 
orthorhombic structure, those perovskites with 
relatively larger B** ions and smaller A*’ ions. An 
intermediate group are those compounds labeled 
pseudocubie which, although not cubic, have not 
been definitely identified as to structure type. 
This group separates the cubic from the orthorhombic 
perovskites in figure 2. Superimposed upon the 
cubic group are those ——— ar -- are tetragonal 
(c/a>1) or pseudo-tetragonal ( ferroelectric 
and antiferroelectric, se il o area showing 
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these compounds can be drawn satisfactorily, as 
» 


indicated in figure 2. However, if the structures of 


various solid solutions are discussed in terms of this 
diagram it is found that there are certain dis- 
crepancies in the boundaries of the ferroelectric field. 
It is true that a complete series of solid solutions 
exist between BaTiO, and PbTiO,, all of which are 
ferroelectric-tetragonal ca i}, Hlowever, com- 
plete solid solution also exists between BaTiO and 
SseTiO . only in this case the solid solutions revert to 
cubic symmetry (at room temperature) with a 
relatively small content of SrTiO,. These properties 
are in contradiction to the ferroelectric field shown 
in figure 2. The properties of the solid solutions 
between PbTiO... PbZrO... and PbHfO,. can be 
explained by introducing a rhombohedral field 
between the tetragonal fields, figure 2. 

The discrepancies ean be partially explained by 
utilizing the theory of polarizability of the ions 
Ignoring the polarizability values for the tetravalent 
ions because not enough values are available to make 
a decent analysis of the results, a three-dimensional 
plot can be made with polarizability of the divalent 
ions as the third coordinate. <A partial cross section 
of this type of diagram is shown in figure 3. In this 
diagram the polarizability of the divalent ions is 
plotted against the radius of the tetravalent ions. 
The diagram correlates the solid solutions of Ba. 
Pb, and Sr titanates. The rhombohedral field has 
been so drawn in figures 2 and 3 to indicate the 
ferroelectric phase observed by Shirane and Hoshino 
(32, 25] in solid solutions of PbZrO, and BaZrQs. 

The symmetries of the ferroelectric solid solutions 
can best be seen if a three-dimensional model is 
constructed with polarizability as the third dimen- 
sion. A diagrammatic representation of such a 
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three-dimensional model is shown in figure 4. Thy 
field of ferroelectricity and antiferroelectricity then 
becomes a volume instead of an area, and _ the 
symmetry of solid solutions is correlated by con- 
necting a straight line between the end-member 
compounds plotted in three dimensions. This 
straight-line relationship probably only holds true 
for solid-solution series that do not cross the cubie- 
orthorhombic boundary. Other solid-solution series 
probably can only be shown by i curved line connect- 
ing the end-member compounds. It is recognized 
that the polarizability of the B™ ion also plays an 
important part in determining the symmetry of a 
perovskite compound. Theoretically the third di- 
mension should be some weighted average of the 
polarizability of both the A** and B** ions. How- 
ever, as the exact nature of such an average cannot 
be calculated at the present time only the values for 
the divalent ions are used. In figure 2 the area of 
ferroelectricity and antiferroelectricity, represented 
by the field of tetragonal ca 1), rhombohedral 
and pseudotetragonal ca 1), can now be con- 
sidered Lo be a two-dimensional projection onto a 
basal plane of a volume of such structures formed 
by plotting polarizability as the third coordinate. 

It should be possible to correlate the solid solutions 
of other perovskite compounds of this group on the 
basis of the proposed diagrams. However, tt 
should be remembered that the boundaries shown 
are only approximate; no account is taken of the 
polarizability of the B*t ion, and the compounds 
distant from each other in ionic radii, or whose 
solid-solution series cross svmmetry boundaries, do 
nol necessarily have a linear relationship 

In the system BaTiO,-CaTiO; [33], the tetragonal 
solid solution extends further than would be shown 
by a straight line connecting the two compounds 
figures 2,3, and 4. The reason for this is either that 
the BaTiO, position is not shown high enough in the 
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field of ferroelectricity because no account is taken 
of the polarizability of the Ti* ion, or that the solid- 
solution S\ mmetries of additions of CaTiO to BaTiO, 
cannot be shown on a straight-line basis because of 
the large difference in ionic radii and polarizability. 
Probably both effects are applicable in this case. , 

The rhombohedral area has been indicated in 
figures 2 and 3 because such a symmetry has actually 
been found in this study. However, the orthorhom- 
bic field suggested by McQuarrie and Behnke [33] 
for the solid-solution series BaTiO,-BaZrO, and 
BaTiO.-CaZrO, has not been included, as the X-ray 
powder patterns of this orthorhombic phase can 
only be indexed as pseudocubic and the symmetry 
is inferred from dielectric data. This orthorhombic 
field seems very logical and probably should be in- 
cluded somewhere in the diagrams. 

McQuarrie [34] has reported a tetragonal phase 
oeeurring in the CaTiO,-SrTiO, system, between 55 
mole percent SrTiO, and 85 mole percent SrTiOs. 
This tetragonal phase is not shown in figures 2 and 
t, although it falls in the pseudocubiec area separating 
the eubic SrTiO,-type structure from the ortho- 
thombic CaTiO, type. These diagrams predict 
that, if this system forms a complete series of solid 
vlutions, there should be a field of symmetry other 
than cubic and orthorhombic. X-ray patterns of 
specimens prepared by McQuarrie, as well as several 
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prepared in this study, do indeed seem to indicate 
tetragonal symmetry. However, as SrZrO;, which 
is also in this pseudocubic area, is definitely not 
tetragonal, some doubt still exists as to the true 
symmetry of this portion of the diagram. It is quite 
likely that all of the area covered by the pseudocubic 
field in the present diagram does not have the same 
symmetry. Perhaps the area grades from tetragonal 
through orthorhombic into still other symmetries. 
It is apparent that much more work remains to be 
done on the symmetry of solid solutions of this type 
before a finished diagram can be presented. 

Theoretically figures 2 and 4, especially the three- 
dimensional diagram (fig. 4), should be usable to 
predict the symmetries of the solid solutions of any 
two or more compounds. However, it must be 
remembered that the positions of the boundaries, as 
shown, are only approximate because of lack of ex- 
perimental data on solid solutions. Therefore, at 
present, the use of figure 4 is limited to selecting the 
most probable of several possible symmetries. In 
many cases, several possibilities may seem to be 
equally probable and experimental evidence will be 
required before a final decision can be made. The 
diagram should prove of considerable help in inter- 
preting inconclusive data that might otherwise 
receive an incorrect or indefinite interpretation. 

The present diagrams, figures 2, 3, and 4, have 
been constructed from the data observed in ceramic 
bodies at room temperature. The appearance of 
these diagrams would be quite different if constructed 
for other temperatures. For instance, at about 
500° C the diagrams would show no field of ferro- 
electric-antiferroelectric compounds and solid solu- 
tions. Between —30° C and 0° C BaTiO, would 
be orthorhombic instead of tetragonal, and the ortho- 
rhombic field, which has been left out in the present 
diagrams, would play a much more important role 
in the diagrams. However, at still lower tempera- 
tures the rhombohedral symmetry of BaTiO; would 
probably increase the area of rhombohedral solid 
solutions at the expense of the orthorhombic field. 
It seems certain, however, that if enough information 
on various compounds and solid solutions were availa- 
ble, a diagram of this sort could be prepared for any 
given temperature, which would explain the solid 
solutions of the perovskite compounds at that tem- 
perature. 


4.3. Mixed Oxides of Trivalent Ions 
a. General 


Mixed oxides of the trivalent ions are known to 
form perovskite-type structures [26, 19]. Very 
recently a patent has been granted B. T. Matthias 
on an electrical device embodying ferroelectric 
lanthanum-containing substances. The compound 
LaAlO,. (and LaGaQO,) was claimed to have ferro- 
electric properties, although no proof of such ferro- 
electricity or any description of the electrical proper- 
ties was given [35]. The question as to whether or 
not these rare-earth aluminates are ferroelectrics is 
very important. In agreement with single-crystal 


results of Matthias [35], rather low dielectric con- 








stants at room temperature were found for ceramic 


specimens of relatively high porosity 


The X-ray patterns of these rare-earth aluminates 
have rhombohedral symmetry, as shown in table 3. 
The composition, structure type, and symmetry, 
where known, for the mixed oxides of trivalent ions 
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are listed in table : 
given only for materials studied in the present work 
and results of previous investigations are listed on} 
when controversial subjects are involved or ¢h, 
particular composition has not been studied in the 
present work. 
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inkuchen, 
ilso observed in present work for powdered material; however, intact 


, the compound BiAlO 


Continued 


, 
References and discussion 

Present work. a=5.545, b=7.851, c=5.562 A. 

Naray-Szabo [26] 

Keith and Roy [19]. “Very slight distortion, 
possibly monoclinic.’’ Orthorhombic. 

W. L. Roth [11], Yakel [43]. Unit cell twice 
that of simple perovskite. 

Keith and Roy [19]. Orthorhombic 

Keith and Roy [19]. Probably orthorhombic 
perovskite. 

Keith and Roy [19]. Probably orthorhombic 
perovskite. 

Geller [54]. Same space group listed as given 
hy Megaw [15] for CaTiOs. 
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Largest distortion from cubic symmetry of 
any of the pervoskites in the present study 
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b. YAIO 


The compound Y AlO, was assigned to a new strue- 
ture type, YCrO;, by Keith and Roy [19]. How- 
ever, the compound YCrQ, has been called a perov- 
skite by Looby and Katz [36]. The X-ray pattern 
for the YCrO,-type structure in Y AlO, specimens can 
definitely be indexed as a perovskite type, as can be 
seen in table 4, and is similar in all respects to the 
CaTiO,-type structure. This structure has the addi- 
tional advantage of having a much smaller unit cell. 
Actually the distortion from a cube is no greater 
than that found in CaZrO, or CaSnO; and does not 
deserve a new structure-type designation. There- 
fore, the orthorhombic-type structure, originally 
assigned to CaTiO;, may now be considered to be 
the more appropriate structure for the compounds 
referred by Keith and Roy [19] to the YCrO, type. 

YAIO; has been reported to form a garnet (solid 
solution) type structure by Keith and Roy [19] who 
state that the garnet structure is a low-temperature 
form and the perovskite (YCrO,) type structure is the 
high-temperature form. This observation does not 
agree with the evidence found in the present study 
table 4). A 1:1 mixture of Y.O, and AIO 
heated for 1 hr at 1,500° C contained about 50 per- 
cent of the garnet-type structure and 50 percent of 
the perovskite-type structure. When this specimen 
was reheated at 1,550° C for 19 hr the amount of the 
perovskite-type structure had _ greatly 
relative to the garnet tvpe. A second specimen of 
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YQ ;:Al,0; heated at 1,835° C ° for 1 hr contained | 
even less of the perovskite structure. It seems that 
attainment of equilibrium between the garnet and 
perovskite structures in YAIO; is a slow process. 
and one of these structures might be only a mete. 
stable phase and have no true equilibrium position, 
A detailed analysis of the phase equilibria throughout 
the Y,O,-Al,O; system would be necessary to answer 
this problem completely. 


4.4. Classification of ABO, Compounds Containing 
Double Oxides of Trivalent Ions 


As very little is known quantitatively about the 
polarizability values of the trivalent ions in peroy- 
skite compounds, this factor has not been used iy 
the present work to describe the classification of 
these structures. All the compounds studied are 
indicated in figure 5 according to the radii of the 
constituent A and B tons. As a complete diagram 
would list each compound twice, the B ion, for con- 
venience, is always taken as the smaller ion. It can 
be seen from this diagram that all of the compounds 
in the upper left of the diagram, large A** and small 
B ions, form perovskite-tvpe structures. This 
group is followed by a belt of Tl.O,-type structures, 
most or all of which probably represent solid solutions 
rather than true compounds. The lower left corner, 
containing small ions in both the A and B positions, 
shows a field of corundum (or ilmenite) type strue- 
tures, probably solid solutions. The upper-right 
corner, in which both A and B are large ions, contains 
a field of La,O,-tvpe structures, probably solid solu- 
tions rather than compounds. This diagram does 
not differ appreciably from that of Keith and Roy 
19]. However, the perovskite field is enlarged to 
accept all those compounds classified by Keith and 
Rov as of the YCrQO, structure. In addition, a line 
is drawn to separate the field of rhombohedral perov- 
skites from the rest of the perovskite compounds, 
which apparently all have the CaTiO,-tvpe structure 
it should be emphasized here that no compounds of 
the A**BV'O tvpe are known to have a simple cubi 
perovskite-type structure. 

One discrepancy in figure 5 might be found for the 
LaCoO, compound reported by Askhan, Fankuchen, 


and Ward [37 Co is listed by Ahrens [30] as 
having a radius of 0.63 <A, essentially the same as 
that of Cr This would place the compound 


LaCoO, well within the field of orthorhombic perov- 
skites, instead of in the rhombohedral field 
should A possible explanation of this diserep- 
ancy, other than assuming incorrect radii or resorting 
to nonstoichiometry, is the possibility of having large 
polarizability in the lanthanum cobaltate compound 
The compound occurring im the specimen of CeO, 
Til do, mentioned earlier. can be explained by this 
diagram if it is assumed that both Ce** and Ti** have 
been reduced to the trivalent state. Tit* has a radius 
of approximately 0.76 A [30] and the compound 
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perovskites, as shown In figure Dd. The rare-earth 
vanadates reported by Wold and Ward [38] can be 
lassified in this YrTOUPINE as probably orthorhombie., 
as V** is listed by Ahrens [30] as having a radius of 
0.74 A. 

Ferroelectricity in this diagram is represented, 
possibly, by the LaAlO, and LaGaO, compounds. 
Itseems likely that if LaAIO, is ferroelectric, then all 
those rare-earth aluminates shown as being rhombo- 
hedral are also ferroelectric. It is not known whether 
an of the other orthorhombic perovskites like 
LaGaQ. are ferroelectric, but no other reports of 
ferroelectricity have been found for this tvpe of com- 
Assuming, however, that both LaAlO, and 
are ferroelectric [35] and that each has a 
svmmetry type, a very interesting possi- 
bility arises for piezoelectric ceramics. It is known 
that close proximity to a morphotropic boundary 
between two ferroelectric solid-solution phases en- 
hances the electrical properties of a ceramic {7). 
Therefore, as solid solutions of LaAlO, and LaGaO 
should cross a morphotropic boundary, interesting 


pound. 
LaGaQO 
different 


nezoeclectric properties micht be obtained. 


5. Summary 


The structures of the perovskite compounds of the 
type A**B“'O, can be correlated by plotting a three- 
dimensional graph with the ionic radii of A*? and 
Bas two of the coordinates and the polarizability 
of the ions as the third dimension. Although it is 
recognized that the polarizability of both the A and 
B ions plays an important part in determining the 
symmetry, only values for the divalent ions have 
been used. A two-dimensional chart, using only the 
ionic radii, shows that the perovskite compounds can 
be divided into orthorhombic and cubic symmetries 
with a belt of pseudocubic compounds separating the 
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two. Superimposed on the field of cubic compounds 
is an area of ferroelectric structures. A_ three- 
dimensional graph, using the polarizability of the 
divalent ions as the third dimension, shows that this 
area is really a projection onto a basal plane of a 
volume in space encompassing ferroelectric and anti- 
ferroelectric compounds and solid solutions. 

A two-dimensional plot using ionic radii has also 
been constructed for double oxides of the trivalent 
ions. This graph shows perovskite compounds 
bordered by TI1,O,-type structures with corundum 
and La,O, structures at the extreme borders. The 
perovskite compounds can be divided into two types. 
Oxides of the larger rare-earth ions with Al,O, form 
rhombohedral perovskites with alpha slightly greater 
than 90°. The rest of the perovskites, including the 
so-called YCrO,-type compounds, can all be related 
to an orthorhombic CaTiO,;-type structure. No ideal 
cubic perovskites are found for the A**B**Q,-type 
compounds. Ferroelectricity in this graph may be 
represented by the rhombohedral area and possibly 
other compounds adjacent to this field. 


Is 


Thanks are due to L. W. Coughanour and B. 
Jaffe for providing many of the specimens, to 5. 
Marzullo for measuring their physical properties, to 
G. F. Rynders for preparing the specimens contain- 
ing vanadium, and to the staff of the Mineral Prod- 
ucts Division of the National Bureau of Standards 
for many helpful discussions. 
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Thermal Expansion of Some Nickel Alloys 


Peter Hidnert 


Data on the linear thermal expansion of some nickel alloys (manganese nickel, Hoskins 


Alloy 667, Inconel, Evanohm, Monel metal, M-M-M< alloy, Illium alloy, and Waspalloy) for 


various temperature ranges between 20 
of Inconel allovs (probably with 
noted in the expansion at about 700° C 
coefficients of expansion of the alloys are 
coefficients varied from 13.1 to 15.510 
cobalt, and iron to nickel-chromium alloys 


1. Introduction 


Data obtained between 1918 and 1956 on the linear 
thermal expansion of some nickel alloys for various 


ranges between 20° and 1,000° C, are presented. 


2. Alloys Investigated 


The samples of nickel alloys were obtained from 
Wilbur B. Driver Co., Newark, N. J., General 
Electric Co., Schenectady, N. Y., Holabird Quarter- 
master Intermediate Depot, Baltimore, Md., Hoskins 
Manufacturing Co., Detroit, Mich., Manning, Max- 
well, and Moore, Inc., Bridgeport, Conn., National 
Advisory Committee for Aeronautics, Washington, 
D. C., National Bureau of Standards, The Standard 
Calorimeter Co. (succeeded by Burgess-Parr Co., 
Freeport, Ill.), and Universal Cyclops Steel Corp., 
Titusville, Pa. The length of each sample used in 
the determinations of linear thermal expansion was 
300 mm. The cross sections of the samples, their 
chemical compositions, and treatments are given in 
table 1. The chemical compositions of some of the 
samples were determined by the Bureau as noted in 
footnotes of table 1. The compositions of the other 
samples were supplied by the companies that pro- 
duced them. 


3. Apparatus 


The micrometric thermal-expansion apparatus de- 
scribed by Hidnert and Souder {1] ‘ was used for the 
determinations of the linear thermal expansion of 
the nickel alloys. The samples were heated in an oil 
bath or an air furnace. To minimize the effects of 
oxidation or scaling at high temperatures, an ob- 
servation wire was placed in a sharp V-groove cut 
around the sample near each end. Observations at 
each temperature were taken under equilibrium 
conditions. 


4. Results and Discussion 


The observations obtained on heating and cooling 
15 samples to various temperatures were plotted, 
end of this paper, 


Figures in brackets indicate the literature references at the 


and 1,000° C 
a supersaturated content of carbon), irregularities were 
and were ascribed to precipitation of carbide. 
tabulated 

° per deg C., 
0 to 24 percent of chromium) are indicated. 
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are given. During the first heating 


The 
For the range 20° to 300° C, these 
The effeets of additions of copper, 


but the observations for only 5 samples are shown 
(figs. 1, 3, 4). The coefficients of expansion in table 
1 were all computed from the original expansion 
curves. 

The “annealed” samples of Inconel,’ containing 
0.02 to 0.11 percent of carbon, were probably in a 
supersaturated condition with respect to carbon on 
account of their relatively fast cooling from 2,050°F 
(1,121°C). On the first heating of these alloys 
during the thermal-expansion determinations (fig. 1), 
slight irregularities were observed at about 500° and 
700°C in the curves of instantaneous coefficients of 
expansion versus temperature (fig. 2). The irregu- 
larities at about 500°C have been attributed [2a], 
although without conclusive evidence, to the forma- 
tion of an ordered structure based upon the com- 
pound NisCr. 

The author believes that the irregularities in 
Inconel at about 700°C (fig. 2) were caused by pre- 
cipitation of carbide. Such irregularities have pre- 
viously been reported by Hidnert [3] for nickel- 
chromium-iron alloys. The irregularities in the 
latter alloys were observed only during the first heat- 
ing on alloys that had previously been cast or 
quenched from a high temperature. ' 

Figure 2 indicates that during the first heating of 
Inconel, the magnitude of the irregularity at 700°C 
increases with increase in the carbon content. 

No irregularities were observed in the curves of 
Inconel on cooling it from 1,000°C to room tempera- 
ture. The author believes that the expansion curves 
of each sample on repeated slow heating and slow 
cooling will closely follow the curve obtained on the 
first cooling from 1,000°C to room temperature. 
From this it follows that the coefficients of expansion 
on repeated heating and cooling for each sample 
should be the same or nearly the same as those given 
in table 1 on cooling. 

Observations on thermal expansion of Monel 
metal (sample 1861) at various temperatures from 
room temperature to 1,000°C are shown in figure 3. 
The expansion was not reversible on heating and 
cooling the sample during the first test, but the ex- 
pansion was reversible during the second test. 


W. J. O'Sullivan, Jr. [2] reported data on thermal expansion, specific heat, 
thermal conductivity, mechanical properties, electrical resistance, and emissivity 
of Inconel containing 0.02 to 0.11 percent of carbon, 
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Figure 4 shows the observations obtained on ther- 
mal expansion of Waspalloy at various temperatures 
from room temperature to 1,000°C 

The coefficients of expansion of the nickel alloys 
containing 76 to 97 percent of nickel agree with the 
corresponding coefficients of nickel within KX 10~° 
per deg C 3 percent 

A comparison of the coefficients of expansion of 
the samples of Monel metal with those for nickel in- 


dicates that the addition of 27 to 29 percent of 
copper to nickel increases the coefficients of expan- 
sion by a maximum of 1.4X107° per deg C (10 


percent) 

The coefficients of expansion of Illium alloy and 
Waspalloy containing appreciable contents of a 
variety of elements are appreciably less than those 
for nickel for various ranges between 20° and 700° C 

Figure 5 summarizes coefficients of expansion of 
nickel-chromium alloys with and without additions 
of other elements, for the range 20° (or O* (‘) to 100° 
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Figure 2. Instantaneous coefficients of linear expansion of 
Inconel alloys containing 0.02 to 0.11 percent of carbon. 


: e. cooling 


, Heating 


C, from the results in table 1 and other investigations 
[3, 4, The straight line represents the rela- 
tien between the coefficients of expansion and the 
chromium content of the nickel-chromium alloys. 
This line indicates that the coefficient of expansion 
increases slightly with increase in the chromium con- 
tent. The addition of 11 percent of cobalt or 7 per- 
cent of copper with other elements (10 percent) to 
nickel-chromium alloys reduced the coefficients of 
expansion appreciably. The effect of the addition 
of iron with or without other elements is indicated in 
the figure. 

Data on the linear thermal expansion of other 
nickel alloys investigated at the Bureau were pub- 
lished in previous papers [3, 8]. 


oO, 6, 7]. 


The author expresses his gratitude to R. K. Kirby 
for his assistance in the determinations on thermal 
expansion of Inconel alloys, Monel metal (sample 
1861), and Waspalloy. 
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Intermolecular Forces in Air’ 
Abraham S. Friedman 


The intermolecular foree constants for a 6-12 Lennard-Jones potential model have been 
computed for air from second virial coefficients derived from experimental pressure-volume- 


temperature measurements. 


Nitrogen: oxy gen 


interaction second virial coefficients and 


interaction foree constants have also been derived. 


The following symbols are used throughout the 
paper: 

7—Compressibility factor, Pv/RT, dimensionless, 

P Absolute pressure, 


»— Molar volume, cm*/mole 
R Universal gas constant, in units of Py - 


T=Absolute temperature. The triple point of water is 
defined as 7'=273.16°K 

B=Second virial coefficient, em*/mole, 

(=Third virial coefficient, em®/mole 2 

)P=Fourth virial coefficient, em */mole 4 

F=Fifth virial coefficient, em "/mole 4, 


Maximum energy of binding between a pair of molecules 
with a Lennard-Jones pote ntial, 
k=Boltzmann’s constant. 
V,= Mole fraction of the 7th component. A 
h=Lennard-Jones interaction parameter, 7 Nad, in units of 
molar volume, 
¥=Avagadro’s number. 
Collision diameter of a pair 
of Lennard-Jon ‘Ss lnteraction, 
{=Angstrom units em 
é Potential energy of two int racting molecules. 
Separation between centers of two interacting molecules. 


of molecules at zero energy 


Collision diameter of a pair of molecules at maximum 
energy, &,. of Lennard-Jones interaction, 
Compressibility factors for dry air have beer 


determined experimentally and are reported else- 
where {1 By fitting the P-v-T data to a Kamer- 
lingh-Onnes virial equation of state 


za!" _ 


R1 ee. os 


B ( 


the virial coefficients B,C, ), and F were determined 
at Various temperatures in the region between the 
critical and 0° C. These are tabulated in table 1. 
The second virial coefficients of nonpolar mixtures 
can be derived, to a good approximation, from the 
force constants of the pure components (2, 3). If we 
assume that the intermolecular forces of the pure 
constituents of dry air (No, Oy, and A) can be de- 
scribed by a Lennard-Jones model, then we may 
define pseudoforce constants for the air mixture: 


‘ 
€ k mi >: A € k 2) 
irch w up} | he Airborne Equipment Division 
Navy Bureau \ 
F l n brack lic l ! neces at the end of this paper 


411810—h57 


| 
| 


| 
| 
| 
| 


and 
2s ’ _— 
bmix=5 TN) iax=5 TN(DINid)*. (3) 
o 0 i 


If we assume that the air mixture has the com- 
position 78.2 percent of N», 20.8 percent of O., and 
1.0 percent of A (ef. [1]), and the N,, O., and A can 
be represented by a 6-12 Lennard-Jones intermolecu- 
lar potential function with force constants: 

(€/k) x 95.93°, 


(d)x, 3.69 A [3, 4] 


(e/k)o, 118.0°, (d)o, 346A [5] 


and 


119.5°, (d),=3.42 A, [6] 


100.8° and b4;,= 


( €p hk) 
then, from eq (2) and (3), (€/k)at, 
60.82 cm/mole. 

These 6-12 Lennard-Jones pseudoforce constants 
can then be used to compute the values of the second 
virial coefficients for air at various temperatures. 
The values thus computed are in good agreement 
with the experimentally determined second virial 
coefficients [1]. 

Table 2 compares the calculated and experimental 
values of B. 

The second virial coefficient of a binary gas mix- 
ture can be expressed in terms of the second virial 
coefficients of the constituent gases and an inter- 
action coefficient related to collisions of the unlike 
molecules. 

Baix= > ON NB y. (4) 
i 


By, 


We may then write for air, noting that B,, 


Bu (0.79)*Bx, + (0.21)?Bo, +0.3318 Bip. (5) 
Rearranging: 
Bye = By,:o,= (Bar — 0.6241 By, —0.0441Bo.)3.014 (6) 


Using the P-v—T data of Friedman et al., for B,,, 
[1], of Friedman for By.[3], and the tabulation of 
Hilsenrath et for Bo,|6], the N,:O, interaction 
second virial coefficients as a function of tempera- 
ture may calculated from eq (6). These are 


al 
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TABLE 1, Virial coefficients of air 


r Be ( D | 
A cm) mole cm*/ mol i nol nol 
273. 1 12. 57 0. 7343 X10 0. 1468 , RR X 11 
200) 19. 30 7646 TNF oll 
200 $8. 20 1 46308 ole 34 
175 55. 54 2 5492 044 O1US6 
1s 76. 36 3923 — 1101 1184 
* The estimated standard ce ition of the second il I , B, i 
from 4 to 15 percent of the value of the coefficients nd lle " ms of tl 
I abuse he ft t 


! large fract 


other coefficients are a 


TARLE 2. Second 


B 
/ 
Caleul | Ex 

A wot 
73. | 13.4 
2A) 19 ) 
20) bm se 
17 ; 4 
1M 74 ‘ 


ogden OTUGE? 


nleractior second a 

core flicier 

" Rs 
A 

4 ~~? 

oth 44 

‘ 

‘ “4 


tabulated in table 3. These interaction virial 
coefficients can be fitted to a Lennard-Jones 6:12 
potential model. The intermolecular force constants 
for the interaction between a nitrogen molecule and 


an oxygen molecule thus computed are: 
€o/k = 106.4° 

b= 61.1 cm*/mole 

d 3.64 A, 


The intermolecular potential energy is plotted as a 
function of the separation of the centers of the nitro- 





150 
100 F 
50+ 
2 3 
£ oO} + + + 
-50} 
-100} 
-150 | 
Fie RI l Intermole cula pote niial between a nit ogen 
cule and an oxrygen molec ile sol ad l rhe 
rhe pot l rey curves for oxygen-oxygen 
nitrogen tions are indicated by the dashed lines 


gen from the oxygen molecules in figure 1. Th 
equation for this curve is 


(2) (9) 7 


3.64/r)"}. 8 


or 
125.6[(3.64/) 


The collision diameter at the maximum energy of 


binding between the pair of interacting molecules 


$08 A 


is 7 dy 2. For nitrogen :oxyvgen collisions, 7 

1} A. S. Friedman, D. A. Gregory, R. Lindsay, and D. Whit 
In preparation). 

[2] D. W hite, a. Ss. Friedman, and H i. Johnstor PI 
Rev. 79, 235, U1ll (1950). 

[3] A. S. Friedman Dissertation, Ohio State Univ 1950 

[4] D. White, A. S. Friedman, and H. L. Johnston, PI 
Rev. 79, 235, U10 (1950). 

[51 L. Holborn and J. Otto, Z. Phvsik 10, 367 (1922 

(6) J. Hilsenrath et al., NBS Circular 564 (1955). 
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Matrices of Spin-Orbit Interaction of the Electron 
Configuration d*s 


W. R. Bozman and R. E. Trees 


The matrices of spin-orbit interaction of the d‘s configuration have been calculated. 
These matrices are for use in calculations to aid in locating levels and making electron 
configuration assignments in Tetand Tar. The calculations were made by using methods 
and tables given by G. Raecah and by the use of tables of the W-coefficients prepared by 
L. C. Biedenharn, and by Obi and coworkers. The matrices have been checked by caleu- 
lating their eigenvalues with the National Bureau of Standards Electronic Computer 
(SEAC 


The problem of locating and classifying the energy levels of a particular atomic system is 
aided considerably if approximate values of the predicted levels can be determined. These 
approximate values are obtained by comparison with corresponding levels in analogous spectra 
and from calculations of the wave functions and energy levels by well-established perturbation 
methods [1]! The first step is to set up the matrices of the important interactions; this is 
done most easily by using methods developed by Racah [2,3]. The elements of the matrices 
are linear combinations of a limited number of radial integrals, the most important ones 
being the Slater integrals, ?*, for the electrostatic interaction and the spin-orbit integrals ¢;. 
The second step, which is outside the scope of this paper, is to assign numerical values to these 
integrals, and then calculate the eigenvalues and eigenvectors of the matrices, which are the 
energy levels and wave functions of the atomic system, respectively. 

Until recently, most calculations were carried out in an LS-coupling approximation, 
neglecting the spin-orbit interaction [1,2,4]. This was mainly because inclusion of the spin- 
orbit interaction greatly increases the size of the matrices, and hand methods of calculating 
the eigenvalues and eigenvectors are ruled out for all but the simpler systems. However, 
electronic digital computers are becoming more available for carrying out these calculations[5,6], 
so that there is now a use for the matrices of spin-orbit interaction. 

These matrices are especially useful at the present time because theoretical studies of 
atomic spectra are mainly directed toward the spectra of the heavy elements and the rare 
earths, in which it would be a poor approximation to neglect spin-orbit interaction [5,6,7,8]. 
This theoretical interest follows the experimental interest in these spectra (in various 
laboratories throughout the world), which is largely directed toward obtaining material for 
inclusion in the National Bureau of Standards series of atomic energy levels publications [9]. 
The present work was carried out to aid the experimental analyses of the first spectra of tech- 
netium and tantalum, currently being carried out for inclusion in the third volume of the 
publication. 

In table 1, the matrices associated with the spin-orbit parameter ¢, are set up for the 
d*s configuration. The rows and columns of the matrices are specified by the name of the 
LS-coupling state, using the usual spectroscopic notation. In addition, a prefixed subscript 
is added to the designation of the parent term ind‘. This is the seniority number and is used 
to distinguish states of the same S and L values in this configuration (the concept of seniority 
number, and this notation, is explained in reference 3). The methods of Racah were used in 
calculating these matrices [2, 3] (the explicit formula has been given by Trees in footnote 9 of 
reference 4). Tables of the W-coefficients, needed for this calculation, have been prepared 
by Biedenharn [10] and by Obi and his coworkers [11]. The matrices have been checked by 
calculating their eigenvalues with the National Bureau of Standards digital computer (SEAC), 
and comparing them with the correct values known from jj-coupling theory (the latter theory 


is outlined in chapter 10 of reference 1). 
wkets indicate the literature reference n page 96. 


95 








10) L. 
fll] Obi, Ishidzu, Horie, Yanagawa, Tanabe, 


[1] E. U. Condon and G. H. 
1951). 

[2] G. Raecah, Phys. Rev. 62, 438 (1942). 

[3] G. Racah, Phys. Rev. 63, 367 (1943 

[4] R. E. Trees, Phys. Rev. 83, 756 (1951 

[5] G. Racah, Lund Lecture, Proc. Rydberg Centennial Conference on Atomic Spectroscopy, Kungl. F ys. Sallskap. 
Handl. N. F. 65, 31 (1955). 

6] R. E. Trees, W. F. Cahill, P. Rabinowitz, J. 

[7] M. Gehatiah, Phys. Rev. 94, 618 (1954 

[8] D. Sperber, Bul. Research Council Israel THE, 437 (1954 

9] C. E. Moore, Atomic energy levels, Cir. 467, National Bureau of Standards, vol. I (1949), vol. IT (1952), 

C. Biedenharn, Tables of the Racah coefficients, Oak Ridge Nat. Lab., ORNL-—1098 (1952). 

and Sato, Tables of the Racah coefficients, Annals of the Tokyo 


Second Series THI, No. 3 (1953); IV, No. 1 (1954 


Shortley, The theory of atomic spectra (Cambridge University Press, London 


Research NBS 55, 335 (1955) RP2639. 


Astronomical Observatory, 


TABI E l. Vat ces of spin-o hit interaction of the d4s configuration 
Che unit of the element 
J=1/2 

w\ 2 is ~ P P P ‘p 4 P ip ‘Pp ‘Dp ‘T) *pD ‘T) *p D 

s 0) 0) \ 4 yb i) 0 ) ) 0 

14 2y 2 } 
Ss 0) 0 _ _& : 0) 0 0) 
} 3 5 3 
7 2 2 ty 14 2y i 2, 10 
P V3 y Y + ’ + V 0) + ’ 0 
3 9 Is 9 rt) 3 

-_ y l4 V2 y 2y7 y 14 2,5 2, 30 
! V0 _ + t 0) f — | — 
3 IS IS ) 4) LD a 

P 0 2v2 : ly 14  2ye | y2 V3 y 35 0 

3 y y y ’ 2 if 

‘ip 0 t 2y i dy lA v2 d yb y 70 y 105 
4 4) 4 4) y | 60 10 
) 6 | 105 70) 

‘T) () 0 () 0 V M + ’ 3 
2 | ! 20 20 

‘'D 0 0 2) 10 2y5 V 35 yO y 105 27 v6 
3 15 () 60 20 20 20) 

2y 30 105 70 6 ‘ 

D 0 0 0 ind () V p05 ’ ; 
5 10 20 20) 0 


96 





20 


0 
1) 


V6 


%) 


= 3/2 
iDD  @D)D  @D)'D 
vid 0 0 
») 
2, 70 
= 0 0 
30 
l V6 V6 
o 2 t 
5 y 50 Vy 50 
0) 20 5 
9, 3 
() oy oy 
3 ) 
0 yo yo 
b tb 
yo | | 
) b 12 
yo | l 
‘) 12 6 
0 y 105 y 105 
12 ) 
0 0 ah A 
20 
yo V5 5 
Ld } ) 
by 50 5 5 
15 15 15 


97 


(D)*D 


y 70 


12 


y 14 


y 105 


GD)*D 


2y 21 


7V6 


20 


7y5 
20 





al 














Z rl 
the IZA 
> oa 
GA Ac 
a) t) 
t Cl \ 
a) t) 
cA ¢Ac 
a) t) 
t C1 AZ 
0) 
() () 
) 
0s \ ‘) 
oY Peet] 
ride O1Z Zz 
( () 
Ol CI Az 
i 
ny. (0) 
t) 
ei () 
t/*G 
0 0) 
0) ) 
1c Fags | 





Q Ag 


oy’ 











tise 


O1Z sz 


Ol 
cols 


OS | 
col’ 


(y6 








Cc | 


CAzg 


Cc] \z 


Oe: 


oy 


yt 


Ct 


Or +s 


OF 
ous 


Cl 
COL AZ 


) i) 
Ct 
| Ag ) 
t) 
ae 0 
Ct 
G1 0 
t) 
9 0 
EAZ 
0) 0) 
OF 
0) QAy 
09 
() Ay 
( 02 
aE AZ OLE 
6 Ol 
eM g/t 
Ot: 
0 ou’ 
Cl 
0 
Cl A 
cl § 
CL AZ | 
OF ¢: 
iran ris 


aed dv(di 


ye (OD$) 
y(t 

Fae 

Ly ol 

[el At 

12(Ai 
Co(ds 
ady(at 
aricer: 
d:(dat 
d:(d} 
d:\d 
dy(dt 
s\d 


98 











= 


ar 









N + 
>= 


0 0 Hy»( HP) 


0 0 (DD 
0 0 (Dd 
0) 0 D.(D1 
0) () ry<(t 











= x Ly( At 
cA iz 1 bs 
$ IS (at 
ci; LA —e 
C IZ 
GA IZ Ae Le 
< 1Z 
‘ G | 
cj At p Ac ‘ 
1% ; 
() a Faga| 
COLAZ 
Ol Ol , 
| od 
é cols 
Ol ()% ; 
' t (154d 
( col \ , 
4) ' 
‘erace! 
| 
ady(at re racer: 
0 0 ry, (PF) 





99 


GD)§D F)*] FE) 4] 1G)2G 1G??G IGG IG)IG SH)? H ‘H)'H 
>D)*D | y2 V2 0 0 0 0 0 () 
| Vil 5 \ 55 
? 
F)*} \e 9 9 ; 2 ; 6 : . 
| | 5 vil V5 \ 55 
i » + ) ) 
iF)" - 2 3 6 (5 12 
| 5 QV 5 55 0) Dy 55 
GyYG 0 : 0 0 — ae bm Vor 
4 15 1d 1d 15 
( ( 0 vil - - 0) () | 55 | Ilys ay 50 ty 5 vel 
6 15 () 15 15 Alb 
Ene 
oft ¢ 
Cr)°Gy 0 yo V5 29 | v 55 = vil <V 6 <\! tne 
{ 15 15 5 10 ; 15 de 
re 
rn] 
55 5 55 | 5 2 ) lf 
OG 0 yoo \ voS \ \ v wal 
( Z i ») it) z( lo 5 
tvp 
ma 
() () O6 ) { 
HH HH ‘) {) ‘) \ . \ pi. v \ a 
lo ) lo ) ) r 
Lal 
\ 
Va 
4 55 5 4 | ' irl 
HH 9 2 n \: ) Vo 2y | . \ 
i” | is i () i) Ol 
| 
D1 
cl 
DOs 
ist 
J 2 
" plit 
ype 
G)'*G H)?H H)*H 1)7] red 
] 139 (dig 
wit 
(a (; ; \ 1O y 26 1) nt 
a ) > H)*H 1)*] m 
qu 
. , Vy 10 y 65 y 50 , | V6 
(H)?H : H)*H do 
) ’ ») b 2 2 
lak 
liv 
2t 6 | 1d : 
oH) = = ; i 1)71 vs 0 
0 i 6 2 
. me 
jv ‘ 
‘T)?] 0 \ pvio 0 pre 
ob b 
Tal 
cos 
WASHINGTON, September 12, 1956. | 


100 





1)*H 


Vol. 58, No. 2, February 1957 Research Paper 2740 


A High-Speed Computer for Predicting 
Radioactive Fallout 
J. H. Wright, L. Taback, and H. K. Skramstad 


A very high-speed analog computer has been developed for predicting the geographic 


pattern of radioactive fallout that will oceur following a nuclear-weapon explosion. 


Input 


data consist of wind magnitudes and directions as functions of initial height in the cloud, 
and the distribution of radioactivity in the cloud as a function of the particle size (in terms 


of settling rate 


and initial height in the cloud. 


Output data are presented on a cathode- 


ray tube, so that the luminance at any point on the sereen represents the intensity of fallout 


accumulated at that point 


1. Introduction 


The computer deseribed in’ this paper was de- 
veloped for the U. S. Weather Bureau and_ the 
\Ibuquerque Operations Office of the Atomic 
Energy Commission, with the immediate objective 
of operation in the Pacific Proving Grounds during 
the 1956 tests of nucleal weapons. Two essentially 
dentical units were built and used by the fallout- 
prediction group as aids for safety purposes. These 
nits, as reported here, were experimental models 
A quite different 
[vpe ol special-purpose analog computer, developed 
bv the Sandia Corp., was also emploved Another 
machine, using mechanical-optical techniques, Is 
inder development at the Los Alamos Scientific 

The most distinctive features of the 
Sureau of Standards machine are. its 
instantaneous speed of solution, straight- 


lacking engineering refinements. 


Laborator\ 
Nationa! 
virtually 
forward procedure for setting in data, and adapt- 
ability to variations of the “‘model’”’ of the radio- 
Various limitations, as well as some 
possible refinements, are discussed. 

In the past, hand methods of solution have been 
ised extensively, These require considerable sim- 
plification of the problem and exploitation of the 
yperator’s judgment to achieve the solution in a 
an hour or less. Very large 
digital computing facilities have also been used 
with considerable reduction in’ the computation 
lime. Both of these techniques will be of continuing 
Importance; the hand methods, because they re- 
quire negligible outlay for equipment and installa- 
tion; the large digital computers, because they can 
do increasingly complicated problems simply by 
taking more time, In contrast, the desired objec- 
lives for the special-purpose computer were: 


active cloud 


reasonable time, e ia 


1. Small size and reasonable portability. 

2. Low cost—under $10,000. 

3. Problems set up without the need for mathe- 
matical training in numerical analysis or digital 
programing. Actually, the operator need not be 
lamiliar even with mathematical functions such as 
cosine or exponential. 

+. Very rapid solution of the problem. 

5. Direct visual presentation of the final solution. 

6. Rapid problem modification. 
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These requirements have generally been met in 
the first model, with some reservations. The solu- 
tion time is considered to be “instantaneous” be- 
cause it is negligible compared to the time required 
for human observation of the output display of 
data. The estimate of cost is crude and is based on 
construction of 3 or 4 duplicate machines, exclusive 
of engineering time or test or of commercial profit, 
The over-all accuracy of the machine has not been 
studied systematically, because of the immediate 
need for field service, and cannot be discussed. The 
serious source of error was the conventional 
tube emploved, which, however, pro- 
vided the most appealing and useful feature of its 
operation in the first application. 

It should be noted at the outset that this computer 
is not a direct simulation type of analog machine, 
but rather a scheme for using analog techniques to 
mechanize the problem in a manner somewhat 
similar to the organization of the problem in a large 
digital computer. In terms of apparatus, this 
machine consists primarily of function generators 
and scanning voltages. ‘Time, as represented in the 
flow of machine operations, has no relation to the 
time variable in the original physical model, but is 
used for sequencing purposes only. 


2. The Fallout Problem 


The assumptions made in the solution of the 
problem are not unusual in the field of fallout cal- 
culations, and will be listed here in some detail. 
First, it is assumed that a cloud of radioactive 
particles has formed, with a definite shape and size, 
and, at some instant considered as zero time for the 
calculations, the various kinds of particles at various 
locations in the cloud begin to settle, being displaced 
laterally by the winds as they fall to the ground. 
This ignores the possibility that some particles, 
notably very large ones, may have reached their 
maximum height and started falling while other 
types of dust particles are still rising toward their 
“initial” positions. The cloud is assumed to have a 
vertical axis of symmetry, whereas actually the 
horizontal winds can tilt and otherwise distort the 
cloud while it is forming. Vertical winds are neg- 
lected. The large-scale vortex motions caused by 
the forces developing the cloud, and any effects on 


most 


cathode-ray 








the wind pattern resulting from the energy of the 
burst, are also neglected. It is assumed that the 
characteristics of the particles are sufficiently 
described by specifying a rate of fall as a function of 
altitude, and a particular weighted value of radio- 
activity for each of the various ranges of fall rate 
and initial height. We will hereafter refer to a 
“slowness” variable, s, because the corresponding 
variable in the machine is proportional to fall time 
rather than rate. The slowness varies greatly with 
particle size. However, it is assumed that the 
relative times spent in the various altitude layers are 
the same for all particles starting from the same 
height. Thus, the law of fall can be made correct 
for some intermediate size of particle, but becomes 
increasingly incorrect toward the extremes of particle 
size. 

Some of the above assumptions become 
ingly important for calculations at positions nearing 
ground zero, the site of the explosion. Eventually, 
perhaps within an area of two mushroom diameters 


increas- 


centered on ground zero, these factors pretty well 
invalidate the detailed computations. The mush- 
room diameter refers to the mean diameter of the 


widest portion of the particular cloud considered. 
A basic assumption is that the wind at each level 
is invariant over the entire area of fallout, and does 
not change during the entire time of fallout. Time 
and space variability of winds therefore introduce 
increasingly large errors for particles taking a long 
time to fall, and for fallout far from ground zero. 
The relative importance of such variations during 
actual operations is the subject of present research 
by Sherman,' Nagler,’ and others. Detailed caleu- 
lations of the effects of variability at all altitudes 
and areas, and during the time of fallout, could 
complicate the problem enormously) by several 
orders of magnitude. Thus, practical considerations 
of operating conditions become all important, as to 
whether the requisite detailed weather data can 
actually be obtained, and soon enough to be valid. 
For example, if a fallout prediction is to be made 12 
hr before detonation of a bomb, and must allow for 
particles requiring 24 hr to reach the ground, a very 
detailed 36-hr weather forecast is such 
forecasting is difficult, if not impossible. Roughly 
speaking, we may say that bevond about a 200-mile 
radius, or 12-hr falltime, fallout predictions neglect- 
ing variability are of doubtful validity, at least 
without extensive reappraisal by weather experts. 
The radial distribution of radioac tivity at any one 
initial height in the cloud is assumed to be the same 
for all types of particles. As will be seen, the present 
computer does not control this radial distribution 
closely. Thus, it may be roughly uniform, or fall 
off with increasing radius, depending on whether or 
not a compensating waveform is applied to the in- 


necessary 


tensity grid of the cathode-ray tube. (An arbitrary 
profile could be generated, but such a refinement 
was considered unimportant in the first model.) 


is assumed that no diffusion occurs during 


Finally, it 


! Los Alamos Scientific Laboratories 
?U.S8. Weather Bureau. 


fallout, so that particles of a given size originating 
in any specified layer of the cloud maintain their 
relative positions as they fall. Thus, the cale ‘ulation 
2 wind effects on partic ‘les on the axis is sufficient: 
i. e., the entire “slice” of cloud (considering partic les 
“ one size range) moves as a piece during the fallout. 

The machine does not need to distinguish between 
the number of particles per unit volume, p;, and the 
activity per particle, qi, for particles of a given size 
(or The machine input known as 
activity, a, is the weighted activity for that size of 
particle for an entire layer of cloud. Thus, a, 
dieiV;, Where the subscript j refers to the jth layer 
of cloud. V, is the volume of this laver, so that 
p:V, is the total number of particles in the jth layer 
of cloud. For a given type of nuclear we apon, values 
must be assigned to the various q;;, p;, and V;, so as to 
define a function of two variables, a=/(h,s), where | 
represents height in the initial cloud, and s the slow- 
ness (time to fall a unit height) typical of a given class 
of particle. This function @ is commonly called the 
“eloud activity model,” In practice, the cloud 
activity model can be supplied to the machine 
operator as a family of curves, and he need not even 
be familiar with mathematical functions. The 
activity function is set up so that the output is in 
terms of percentage of total activity in the cloud, 
The family of curves are adjusted with A as a parame- 
ter for 0, 20, 40. 60. 80. and 100 percent ° of the 
maximum altitude of the cloud. This normalizing 
procedure provides for solutions over a considerable 
range of weapon sizes without readjustment of the 
cloud activity model as long as the shapes of the 
clouds remain the same. 

It must be noted that the machine does not make 
any detailed allowance for decay in radioactivity, 
either during the time of fall afterward. The 
assignment of a;; to a given range of particles at a 
given initial height must include the allowance for 
decay. Furthermore, for a particular activity model 
set up in the machine, the operator cannot make a 
simple single adjustment to show the fallout at some 
arbitrary earlier stage, 1. before all the particles 
included in the model have reached the ground. 
However, the activity model caa be modified to 
include only those particles that will reach the ground 
before the specified time. In any the model 
can be set up on either of the following bases: 

Integrated dose. The total from time of 
deposition to some specified time, such as 48 hr., ¢ 
10 years, after the burst. Often the “infinite dose,” 
for complete decay, is used in estimating possibility 
of hazard to permanent residents in an area. The 
unit of measure is the roentgen, R, or milliroentgen, 
mr. 

Dose rate. The intensity of radiation at some 
selected time after fallout, such as 48 hr., This 
quantity may be expressed in roentgens or mill- 
roentgens per hour, and is appropriate to considera- 
tion of evacuating personnel or exploring a con- 
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It is a little more convenient than 


taminated area. 
the integrated dose in deriving a model, because the 


time of arrival need not be known. In practice, it is 
common to use one law of decay applying to all the 
batches of particles, the gross fission-product decay 
Agt™'*, in which A is the dose rate 
time ¢ after the explosion. Ap is the 
dose rate at + ni time, such as roentgens per hour at 
1 hr. On this basis, the completed fallout celeula- 
tions for dose rates, at some time such as 48 hr. after 
the explosion, can be converted by a single factor to 
provide the dose rates at some other time, such as 
10 days later. 

Despite this simplification of the problem, the 
required amount of computation for a nearly in- 


e xpression, - 
measured ¢ 


stantaneous presentation is quite formidable. To 
illustrate this, let us consider the following break- 


down of the problem into discrete elements appro- 
priate to digital computation. The cloud is divided 
into 20 layers vertically, and 25 elementary 
sectional areas in each layer. o> there are 500 
elemental volumes in the cloud, each containing a 
distribution of particle SIZES, with attendant varia- 
tions in “slowness” of fall. a we divide this range 
of particles into 10 classes, there are effectively 5,000 
lumps of material in the deat. each requiring calcu- 
lation of displacement by the wind in each wind 
laver through which it falls. On the average, 10 
multiplications must be performed to obtain the final 
horizontal displacement of each elementary lump of 
material as it falls to the ground, each the product 
of a time of fall and a wind vector. If the solution 


crTross- 


is to be completed within 1/10 sec, the required com- 
puting rates are: 

(a) 50.000 function table values per second, for 
radioactivity. 

(b) 500,000 simple multiplications per second, for 


partial displacements. 
Appropriate combination of 
results, and presentation on a map. 

The display device would be required to generate 
“picture elements” of varying intensity at the rate 
of 50,000/sec, one for each representative lump of 
radioactive dust at its final fallout location. These 
picture elements tend to be superimposed, because 
light ery les originating in a given laver will fall 
on top of heavier particles originating from lower 
altitudes in positions relatively “downwind” with 
respect to the light particles. This situation will be 
referred to as overlap. Thus, the display device 
must be able to accumlate these overlapping contribu- 
tions, if the computer proper is to be relieved of this 
burden. 

The accuracy requirements on individual opera- 
tions are rather low, generally ranging from errors of 
2to 10 percent. An analog technique that endeavors 
to meet approximately the requirements of the 
example will now be described. 


(¢) these partial 


3. Wind Computations 


_ Let us assume that the mushroom cloud (see fig. 1) 
is divided into horizontal lavers at graded intervals of 





Figure 1. Mushroom cloud. 





Wind hodograph. 


FIGURE 2. 


altitude such that a particle of typical size spends 


equal times falling through successive layers. Let 
h represent the initial height of any particle. Now 


let us define a variable h*, a function of h, which is 
proportional to the time for a particle to fall from 
any specified height to the ground. The variable 
h* has the value 1 at the top of the bottom layer, and 
n at the top of the nth layer. Let V; represent the 
average speed and 6, the average wind direction in 
the 7th layer, @ being measured clockwise from north. 
The cloud diameter as a function of height is repre- 
sented by d(h) 

Figure 2 shows a wind hodograph. The dotted 
line may be considered to represent the horizontal 
projection of the path of a hypothetical weather 
balloon having a constant rate of ascent. If an 
appropriate altitude scale is marked along such a 
curve, a vector such as V, joining altitude points 
h, and h, would represent the direction and magnitude 
of the average horizontal wind exerted on this 
hypothetical balloon as it rises completely through 
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the laver between heights h, and hy. Such hodo- 
graphs are commonly used to obtain measurements 
of winds aloft. In the case of the descending motion 
of particles, the hodograph would be drawn in the 
reverse sequence, starting with V, instead of 4,. 
However, the total relative displacement would be 
the same, if the time spent in each laver were the 
same as for the balloon. The horizontal displace- 
ment of a particle while it falls from fA, to A, is 
obtained by multiplying this wind vector by the 
time required for this particular particle to fall 
through this layer. A particle of size @ falling from 
the top of the third laver to the ground would be 
displaced in an east-west direction by the amount 


sin @.], 


r=s(a)[V, sin 6,4 V2 sin 0+ V 
where s(a) is the slowness for this size of particle, 
the time to fall any one height interval. <A similar 
cosine expression holds for the north-south compo- 
nent of displacement. 

It is apparent that there is an approximation 
introduced by assuming a finite number of wind 
lavers, because the dotted line of the wind hodograph 
deviates from the straight-line vectors joining the 
laver boundary points f,, hy, However, this 
introduces no error for the particles falling com- 
pletely through a given laver, because the total effect 
of the wind is specified correctly. There is a small 
error, within the laver only, for those particles that 
start to fall from some height within the particular 
laver under consideration, because the analog com- 
puter interpolates along the straight line in space 
If a digital computer were used, it would be necessary 
to divide the cloud into a great number of lavers, 
or to provide an interpolation routine 

The basic equations for the final ground position 
of a particle of slowness s,, originating on the axis 
of the cloud within a laver for which n<h*<n+1, 


are therefore 
7 | S5¥% sino h*—n)V,., sin 0 | 

1 
cos } 


~— 
y , | av cos 4 hh n io 
where x represents the fall position east of the 


explosion site and y the position north of the explo- 
sion site. In each equation, the second term, such 
as (h*—n V, , sin 6... represents the initial displace- 
ment of the particle while falling through a fraction 
of laver n 1) enroute to the top of laver n. 


etc. 


4. Description of the Computer 


A simple block diagram of the apparatus is shown 
in figure 3. From the viewpoint of the operator, 
there are four basic units: 

(a) The winds unit, with knobs for setting in the 
magnitudes and directions of the mean winds in 


each of the 20 lavers. 
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b) The cloud geometry unit, including the 


cloud diameter and cloud raster functions. There 
are several knobs for adjusting the cloud shape, and 
one knob for setting in an over-all scale factor for 
cloud Size. 

(c) The cloud activity unit, which has 36 knobs 
for setting in the “model” of distribution of radio- 
activity in the cloud. 

dl The display unit, which has J and Y scale-factor 
controls (marked y-amp and s-amp), a brightness- 
level control, and 2 brightness-slope control. 

There is also a central control unit containing all 
timing and circuits, summing amplifiers, 
slowness modulators, and final « and y= driver 
amplifiers. It needs no attention in normal opera- 
tion. In addition, there is a photometer used to 
measure the brightness on the direct 
It is an electrically and mechan- 


sweep 


screeh as a 
measure of fallout 
ically separate entity. 

For explanation, let us first that the 
cloud is collapsed laterally so that all the particles 
are on the vertical axis of the cloud. Furthermore, 
consider that we are dealing only with the /ightest 


consider 


1. €.. Slowest particle S give nmany consideration in the 
probl m., We will later modify the solution to include 
heavier particles, and those off the axis of the cloud 
A periodic saw tooth voltage is venerated to represent 
the variable A*. This is called the height sweep. 
The voltage rises linearly with time to a nominal 

100 v, is reset to zero in about 40 usec, and repeats 
each 1/10 see, seanning through the full range of h* 
each time. This is the basic cvcle of computer 
operation, and the entire solution is generated anew 
with each sweep through the full cloud height. This 
voltage is applied to the winds unit, a function 
generator with 20 elements sequentially activated 
as the input voltage rises, but with their outputs 
added together at all times. At the start, with the 
h* voltage at 0, none of the elements produces an 
output. As the A* voltage rises from 0 to a nominal 
value of 5 v, the output of the No. 1 element, 
representing the lowest wind layer, also rises from 


0 to +5 v. As h* rises from +5 to +10 v, the 
output of element 2 becomes activated, and rises 
from 0 to +5 v. thereafter remaining at +5. This 


goes on until 4* reaches a nominal value of + 100 ¥, 
and all 20 function-generator elements are pro- 
ducing outputs. 
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Frat RE 4. Schematic of winds unit, 


This sequential and cumulative action provides 
the necessary effect of only 1 wind contributing to 
the total displacement of particles falling from the 
top of the first layer, but 20 winds contributing 
to the total effect on particles originating at the 
very top of the cloud. For a particle originating, 
say, halfway in layer n, the cumulative wind effect 
is represented by the sum of the earlier sustained 
voltages, representing lower altitude winds, plus 
half the full scale voltage from the partly activated 
nth element. The situation is shown in figure 4 for 
particles in the first five lavers of altitude. Actually, 
before these voltages are summed, they pass through 
attenuating potentiometers, adjusted according to 
the mean magnitude of the wind in each layer. The 
nominal +5 v represents a maximum-scale wind, 

g., 100 knots. If there is no wind in a given layer, 
that potentiometer is simply zero. The 
voltages are also resolved into + and y components 
by sine-cosine potentiometers set according to the 
wind direction, for each layer. 

Figure also provides a simplified electrical dia- 
gram of one element, for wind laver. When h* 
is 0 at the beginning of the sweep, the lower diode, 
connected to a nominal 5 v (about 4 v is used, 
to allow for forward drop in - diodes), is non- 


set to 


one 


one 


conducting because it is held in a reverse bias con- 
dition by the action of the Soe negative voltage 
V. The uppe r diode is conducting, meanwhile, 


and sets the input to the wind potentiometer at 0 v. 
The anode of this diode is shown connected to +0.6 
Vv, a value adjusted j operation to compensate for 
the diode forward drop so that the wind-potentiom- 
eter output actually is 0. As the A* voltage rises, 
the voltage at the arm of the altitude breakpoint 
potentiometer will at some time, depending on its 
_s above 0 \ While is between 0 and 
both diode ‘sare none ondue ting, and 2 voltage 
buil Ie up across the wind potentiometer, reac hing ‘ 
hominal value of +5 v before being limited by con- 
duction in the lower diode. 
The output of the wind potentiometer drives one 
side of a sine-cosine potentiometer. An intermediate 


rise 


buffer amplifier is necessary to prevent electrical 
loading of the wind potentiometer, and a 1:1 inverter 
amplifier is necessary to drive the other side of the 
sine-cosine potentiometer. AR, is necessary for elec- 
trical reasons to prevent excessive current being 
drawn from the breakpoint potentiometer at high 
settings, and to standardize the forward current 
through the diodes. Resistors R, are identical 
summing resistors feeding ordinary analog summing 
amplifiers, a total of 20 connected to the east-west 
summer, and 20 to the north-south summer. The 
altitude breakpoints are set to occur at 20 steps of 
the h* voltage, equally separated from 0 to +95 v 
(for an A* amplitude of 100 v). The value of the 
nominal +5 v supply is adjusted to make the indi- 
vidual steps add up properly. The 20 rheostats 
labeled R; are calibration adjustments to correct for 
component variations, thus making the full-scale 
contributions from all the wind potentiometers equal. 


As front panel adjustments, the R; rheostats could 
be used to weight the winds so that the 20 layers are 
equally spaced heights, rather than irregular height 
intervals spaced at equal intervals of fall time (for a 
selected particle size). This may prove to be de- 
sirable for convenience and flexibility. If this is 
done, the settings must be changed when shifting to 
a problem with quite different initial cloud height, or 
choice of typical particle. 

Additional sawtooth voltages are generated by the 
particle slowness modulators. Each modulator pro- 
duces four sawteeth while the height sweep scans 
through each cloud layer. The amplitudes of these 
sweeps are scaled according to the relative time of 
fall. At the midpoint of the sweep, for example, the 
vector sum of the outputs of the z and y channels is 
half the full amount. The result is to trace a radial 


line, or “spoke,” that is the locus of fallout of 
particles of various kinds up to the slowest-falling 
tvpe. These spokes are nominally constant-height 


lines, but are curved somewhat, because there is a 
little variation in height, and hence in mean wind, 
during the sean of particle slowness. The effect is to 
multiply each component by the sawtooth voltage, 
but no conventional multipliers are used. The effect 
is achieved by use of voltage-controlled sweep ampli- 
tudes, of whatever polarity is dictated by the wind 
pattern. These modulators are actually integrating 
circuits, using operational amplifiers. They are 
reset SOO times per second by triode reset circuits 
that require about 

To take into account the diameter of the cloud, the 
spot on the cathode-ray tube must be effectively 
enlarged to have a diameter proportional to that of 
the cloud, changing with the voltage A*, as the solu- 
tion of the problem progresses. To spread the spot 
sufficiently, it has been necessary to resort to a high- 
frequency spiral raster, with a circumferential fre- 
quency of about 100 ke and a repetition rate (raster 
rate) that varies from 4 to 20 ke or more. The high- 
est frequency is associated with the smallest spiral 
diameter, because the spiral advances outward at a 
fixed rate until it reaches a specified voltage repre- 
senting the cloud diameter. Thus, larger raster 


20 usec. 
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diameters will be repeated at a lesser rate, and will 
have more turns in the spiral. The computer unit is 
called the cloud raster generator. It is driven by a 
voltage from the cloud diameter generator, which is 
a conventional diode function generator approxi- 
mating a function of one variable (the h* voltage) by 
means of connected line segments having adjustable 
slopes and breakpoints. The z and y voltage com- 
ponents of the spiral are inserted into the final 2 and 
y driver amplifiers as increments to the basic fallout 
positioning voltages. 

The cloud activity generator receives the A* and s 
voltages and continuously generates an output that 
is a function of these two variables. The output of 
this unit can be displayed on the main cathode-ray 
tube when desired, in the form of a family of curves 
of activity @ versus slowness with h* as a parameter. 
Actually, 80 curves appear in the family, but the 
operator sets in curves only at the top of the cloud, 
ground level, and four intermediate points at equal 
intervals of h*. Each curve is made up of five line 
segments. This generator must be unusually fast in 
operation, because s sweeps through its entire range 
of 50 vin 1/800 sec. Furthermore, it should be possi- 
ble to adjust to arbitrary functional variations be- 
tween curves in the family, rather than use poly- 
nomial methods. It is desirable that the operator 
simply adjust a row of knobs until the proper func- 
tion of activity versus slowness is obtained at ground 
level, then independently adjust to match a specified 
curve at the next selected height in the cloud, and so 
on to the top of the cloud. 

Function generators developed by the George A. 
Philbrick Co. were emploved in the first two com- 
puters. They permit such independent arbitrary 
adjustment. However, the interpolation between 
curves for a constant value of the s variable is not 
linear. For example, in progressing from a_hori- 
zontal line to one with a V-shaped notch, the inter- 
mediate curves are truncated, with flat sections, 
rather than V-shapes of intermediate height. The 
device is all electronic, using a matrix of selector 
diodes. Better interpolation is and an 
electronic function generator * recently developed by 
L. Taback of the Bureau may find application in later 
models of this computer. The Taback grenerator 
uses a time-division technique for interpolation be- 
tween curves in a family. 


‘eal 
desired, 


The output voltage of the activity generator is 
connected to the intensity grid of the final display 
cathode-ray tube. It thereby modulates the elec- 
tron beam so as to make the brightness proportional 
to the activity of the particle type at the initial height 
under consideration at any instant. The intensity 
level of the cathode-ray tube not remain 
sufficiently constant (at a specified grid test voltage) 
for a permanent type of calibration. Furthermore, 
the luminance versus grid voltage curve departs 
considerably from linearity in some regions of opera- 


does 


tion. Hence, slope and brightness-level controls 
4‘ L. Taback, A function generator for two inde} ( I 


3, 129 (1956). 


h™, SAWTOOTH VOLTAGE , PERIOD 1/10 sec 
STANDARD AMPLITUDE . 





$s ;,SLOWNESS", ACTUAL PERIOD 1/800 sec 
NOT AS SHOWN ' 


NPUT TO WIND POT. NO.2 
TPUT OF WIND POT NO.2, 2 


3 MAXIMUM 


NORTH OUTPUT OF SIN-COS POT. NO.2 


EAST OUTPUT OF SIN—COS POT. NO.2 


OUTPUT OF EAST WIND SUMMER 


OUTPUT OF EAST PARTICLE SLOWNESS 
MODULATOR 


aeWWw 


Ficure 5. Basic waveforms 





FIGURE 6 


Oulput of north-south winds 


summer, 


are provided. The slope control is an attenuator 
in the activity signal line. These controls must be 
adjusted with the entire computer in operation after 
insertion of a calibration problem with known winds, 
cloud geometry, and activity model. The procedure 
is to first set the « and y deflection factors for the 
cathode-ray display so as to make the display read 
correctly in miles for the eastward and northward 
limits of the pattern. For this purpose the test 
problem can be made very simple. Then the over-all 
size factor for the cloud is adjusted to make the 
mushroom diameter appear to be the proper num- 
ber of miles on the screen. The slope and _ bright- 
ness-level controls for the cathode-ray tube are then 
adjusted to make the photometer read as correctly 
as possible in areas with various ranges of brightness. 

Some of the basic waveforms in the computer are 
shown in figure 5. <All are exhibited on the same 
time scale synchronous with the h* voltage, but to 
simplify the diagram, only 5 wind channels are shown 
instead of the actual 20, and only 10 s sweeps instead 
of the actual 80 during each h* sweep. Four of the 
waveforms are associated only with wind 2, while 
the two waveforms at the bottom show the contribu- 
tions from channels 1, 3, 4, and 5 as well. 

Figures 6 to 14 were taken before the first com- 
puter was completed. They do not exactly repre- 
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FIGURI 7. O itput of no th-south slowness modulator. 


FIGURE 8. 


sent normal operation, providing less than the nor- 
mal 80 for example. Because it was not 
possible to take pictures subsequently, these are 


‘“‘spokes,”’ 


presented to aid the discussion. 

Figure 6 is an oscillograph trace of the output of 
the north-south summer versus h* voltage. Figure 
7 then shows the result after the waveform shown 
n figure 6 is used to modulate the particle-slowness 
voltage. The triangular waveforms actually have an 
envelope that follows the waveform of figure 6. 

Figure 8 is an oscillograph picture of the output 
f the activity for convenience of 
lisplay, and not to represent any actual explosion 
loud. The represent particle slowness, 
and the ordinates represent the output activity, a. 
Each of the individual curves corresponds to a differ- 
nt value of A*. 

Figure 9 shows the display of fallout on the cath- 
de-rav-tube screen in the absence of any intensity 
modulation by the activity generator, and without 
the spreading of the beam by the cloud raster genera- 
tor. Each spoke is the result of one slowness sweep. 
Fallout from higher initial lavers 1 
the cloud corresponds to progression of the spokes 
na clockwise sense for this The number of 
spokes is determined by the ratio of the s sweep 
irequency to the h* sweep Trequency. 


venerator, set 


abscissas 


successively, 


Case. 


Figure 10 was photographed under the same com- 


puter settings but includes the intensity modula- 


tion according to the function a. 


RE 9. Simple spoke display. 


FicurRe 10. Intensity-modulated spokes. 


Figure 11 is similar to figure 10, except for more 
drastic modulation of the intensity. Here only the 
tips of the spokes show, as would be the case if 
only the slowest particles were radioactive and 
numerous enough to be of importance. 

Figure 12 shows the result of superposing the 
cloud raster on the display of figure 11. In actual 
use, the oscilloscope spot is defocused so that the 
spiral pattern does not show clearly. (The appear- 
ance of spirals of widely different diameters on the 
same centers is not proper, and is the result of 
faulty synchronization existing at the time of this 
early picture. 
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Figure 13 shows how the final display appears with | generator. The -I-in. cathode-ray display unit is 
all units of the computer in action The computer mounted on a table at the right On top of it is the 
Was adjusted to produce this apparent fallout with- meter of the photometer, a modified commercial unit 
out knowledge of the characteristics of any actual with a probe area °y in. in diameter The probe IS al 
bomb, hence the pattern is not representative the end of a cable. The cabinet at the left contains 


Figure 14 shows one of the computers. The | power supplies, test equipment, and a spare chassis 
cabinet on the right contains the computer proper, | The second large cabinet is unnecessary. and _ the 


together with an auxiliary 5-in oscilloscope. The 4 | essential] power supplies can be mounted in a small 
chassis below the 5-in scope form the winds unit, unit forming a pedestal for the display unit. Foreed- 
each containing 5 wind channels. Below these. mov. air cooling is employed, with lateral airflow from 
ing downward, are the cloud activity unit, central front to rear over the chassis The total a-e power 
control unit, cloud diameter unit, and cloud raster requirement is about 1.500w. 
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5. Comments on Operation 


Because the machine is used intermittently, and 
accuracy requirements are relatively low, chopper 
stabilization was not employed in the operational 
amplifiers. Thus, balance adjustments are neces- 
sary. ‘These are made after 1 hr. warmup time, the 
rocess taking about 5 min. Following this, some 
readjustment is necessary at intervals of 2 to 6 hr. 
Checks on the brightness level of the scope and 
balance of the activity generator are necessary at 
about the same intervals. A few other basic elec- 
trical adjustments are made at intervals of a few 
davs, including sweep amplitudes and_ altitude 
breakpoints. 

In use, provided the cloud geometry and activity 
model had already been inserted in the machine, a 
completely new pattern of winds could be inserted in 
about 2 min. once winds were specified at the proper 
altitudes. Results of problem modification were im- 
mediately observable on the display, but 2 to 5 min. 
were required to take a useful pattern of measure- 
ments with the photometer probe. 

Modification of the activity model from a set of 
tabular data could be done within 5 min. The shape 
of the cloud could be altered as desired, but it seldom 
needs serious adjustment. The over-all size of the 
doud could be altered instantly by turning a single 
knob. 

Relative position of eround zero on the screen 
could be adjusted as quickly as desired, but with 
some limitations due to shortcomings of the oscillo- 
scope. Likewise, the picture could not be expanded 
greatly, because of the limitation in cathode-ray 
deflection speed to about 0.6 in./usec. This meant 
that the mushroom could not appear as more than 
1% in. in diameter on the screen, due to the periph- 
eral rate of 100,000 turns/sec in the spiral. An 
overlay map on transparent backing was used over 
the display screen. 


6. Cathode-Ray Tube Problems 


Because of “overlap,” any elementary area on the 
eathode-ray screen may receive many contributions 
of beam current during any one computation cycle. 
After one “‘strike’’ by the beam, others will occur as 
the spiral raster moves on out along a given spoke, 
oralong an adjacent spoke. It might seem desirable 
touse a phosphor screen having long persistence, so 
that it would integrate these contributions. How- 
ever, the accumulated light-energy output versus 
eectrical charge resulting from the various strikes is 
not a correct integral. The relation between beam 
current, exposure time, and previous history is quite 
complicated and unsatisfactory. Storage tubes and 
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dark trace tubes were considered, but the practical 
choice was a short-persistence tube. Ideally, the 
persistence should be negligible. The summation is 
accomplished in the averaging action of the pho- 
tometer deflection. Because visual display was 
desired, the computer repetition rate had to be fast 
enough to reduce the flicker in visual effect. A 
compromise of 10 cps was made in this computer, 
although 30 eps would be necessary for really sub- 
stantial reduction of flicker. Actually, the phosphor 
was a fairly severe compromise, being a P-—4 type. 
Theoretically, a considerably faster type would be 
desirable. 


7. Refinements 


Fairly obvious engineering refinements will reduce 
size and improve stability, particularly in choice of 
components. The number of components can be 
reduced, and the winds unit, in particular, could be 
transistorized. Zener voltage diodes can be used to 
simplify the winds unit and eliminate the low- 
voltage precision power supply. 

Greater accuracy and usefulness under some con- 
ditions could be achieved by eliminating the cathode- 
ray tube altogether. The same basic circuits would 
be used, but could be run much slower, thus reducing 
the basic problem of achieving a sufficiently high 
raster rate relative to the slowness sweep. The out- 
put of the activity generator would be gated into an 
averaging circuit when and only when z and y simul- 
taneously have values falling within a_ specified 
elementary geographical area. For any such elec- 
trically probed area, the gate would open many 
times during each solution time, due to the overlap- 
ping previously discussed. The output would be 
read directly on a meter, or could be printed out on a 
standard electric typewriter operating from a rela- 
tively slow analog-digital converter. The presenta- 
tion could have the form of a map, each advance of 
the typewriter representing a standard number of 
miles in a routine scan of the fallout area. 

The present computer design could be extended to 
accommodate more elaborate formulation of the 
problem, such as time or space variability of winds. 
However, it is not known whether it is entirely de- 
sirable to attempt this. 


are indebted to Lester Machta of the 
Bureau, who initiated the project on 
behalf of the Atomic Energy Commission, and to 
K. M. Nagler, also of the U.S. Weather Bureau, for 
his assistance during the project. 


The authors 
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WASHINGTON, September 13, 1956. 
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Research Paper 2741 


Standard lonization-Chamber Requirements for 250- to 


500-Kilovo 


It X-Rays’ 


H. O. Wyckoff and F. S. Kirn ” 


Information is provided on the correction 
dose in roentgens of 250- to 500-kilovolt X-rays 
chamber. Special emphasis is placed on the de 
over-all accuracy of 0.5 percent can be obtaines 


1. Introduction 


In 1953, the International Commission on Radio- 
logical Units recommended |{1] the continued use of 
the roentgen up to 3 Mev for the measurement of 
X- and gamma-ray The physical require- 
ments of standard ionization chambers for measur- 
ing roentgens in the X-ray region up to 250 kv have 
already been reported |2,3,4]. The calibration of 
secondary instruments in at least two laboratories 
also seems to be in agreement [5] for cobalt-60 and 
radium, on the assumption of a fixed value for the 
emission constant (roentgens per milligram hour at 
jem) of radium. However, there is some question 
about the value of the emission constant for radium 
6,7] when its determination is based upon measure- 
ments with small cavity ionization chambers. In 
addition, many of the commercial secondary instru- 
ments show a difference in calibration of the order 
of 20 percent between 250-kv X-rays and 
cobalt, so an interpolation of the calibration for the 
300- to 400-kv X-ray therapy range seems risky. 
Extension of the range of standard instruments, at 
least up to cobalt energies, therefore seems desirable. 
The present report contains the design criteria for 
the construction of a standard chamber for 250- to 
500-kv X-rays. 

A plan view of one type of a standard free-air 
chamber is shown schematically in figure 1. A 
parallel-plate system is contained in a radiation- 
shielded box, B. A known area of the horizontal 
X-ray beam is defined by the diaphragm, D, so 
that a beam of photons passes centrally between the 
plates. A high potential (field strength of the order 
of 100 v/em) on plate H, with respect to the other 
plates, sweeps out the ionization produced in the air 
between the plates. The ionization is measured for 
alength, 1, of beam determined by the limiting lines 
of force to the edges of the collector, C. These lines 
are made straight and perpendicular to the collector 
by the guard plates, G, and auxiliary wires or strips, 
W. The latter are connected a resistance- 
dividing network to grade the potential uniformly 
across the end, top, and bottom gaps between C or 


dose. 
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G, and H. Generally guard plates are not con- 
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factor required for the determination of the 
from measurement with a free-air ionization 
termination of each of the factors so that an 


1. 
sidered necessary above and below the collector 
because there the ionization density is low. 

A measurement in roentgens requires, according to 
the definition [1], that one measure all of the ioniza- 
tion produced in air by the high-speed electrons 
that are themselves produced within the defined 
mass of air. Real difficulties arise in the measure- 
ment of the ionization as required by the definition, 
so the principle of electronic equilibrium is always 
used in practice. 

Within a region under uniform irradiation, ac- 
cording to the principle of electronic equilibrium, 
the ionization produced outside of a specified mass, 
m, by high-speed electrons generated inside m is 
compensated by ionization in m produced by high- 
speed electrons generated outside of m. Thus 
electrons produced out to a distance from the mass 
equal to the electron range contribute to the ioniza- 
tion in the specified mass. Therefore in a standard 
chamber, the guard plates not only must be long 
enough to eliminate field distortion but also long 
enough to assure electronic equilibrium in the collect- 
ing region. For the latter, a distance equal to ralf 
the plate separation has been suggested [4] for X-rays 
up to 250 kv. However, for the higher-energy 
photons used here, the range of the electrons in the 
photon direction may be larger than the sideways 
range. A determination of the distance required for 
electronic equilibrium is therefore desirable. 

Ionization is produced in the collecting region not 
only by the primary electrons, but also by electrons 
produced by secondary (scattered and fluorescent) 
photons. As this secondary ionization is not in- 
cluded in the roentgen definition, a correction is 
required to the measured ionization. However, this 
excess ionization tends to compensate for loss of 
primary ionization resulting from too small a plate 
separation and height, so that the required electrode 
separations are reduced by this effect. 

In the calibration of a radiation instrument by a 
standard chamber, a substitution method is used. 
The exposure dose rate is determined from measure- 
ments of the ionization in the standard chamber, this 
chamber is removed, and the instrument being cali- 
brated is placed at the position formerly occupied 
by the aperture of the standard chamber. There- 
fore, it is necessary to compute the exposure dose 
rate at the aperture from the ionization measure- 
ments. The factors required for such a conversion 
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FIGURE 1. Schematic view of free -(il wonization chamber 
showing a horizontal section through the radiation beam and 
chamber. 
B, Shielded box; C, collector (length, L); D, diaphragm; G, guard plates; H, 


high-voltage plate; W, guard wires or strips; 5, plate separation 


include a knowledge of the weighted position of 
origin of the high-speed electrons relative to the 
collecting volume and the air absorption between 
this position and the aperture. 

Thus the information required for proper 
design of a primary standard chamber (1) the 
plate spacing, S (fig. 1), required to assure dissipation 
of all of the energy of the secondary electrons in the 
air, (2) the criteria for producing a defined electric 
field so that this ionization is collected from a known 
length, L, of the beam, (3) the correction for the 
ionization produced by scattered photons, (4) the 
weighted position of origin of the high-speed electrons 
relative to the collecting region, and (5) the air 
absorption between this origin and the diaphragm. 
Data for item (2) have been obtained by Miller and 
Kennedy for the two-dimensional case [8]. Items 
(1) and (3) have been considered by Attix and 
DeLaVergne [4], by Kustner [2], and by Block [13], 
and item (1) by Ke mp [3], for 50- to 250-kv X-rays. 
The present paper gives information on items (1), 
(3), (4), and (5) for 250- to 500-kv X-rays and 
moderate filtration. Work at higher energies is now 
under way at the Bureau. 


the 
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2. Apparatus 


A vertical, continuously pumped, end-grounded, 
59-section X-ray tube was used. It was connected 
to the ground, center, and top of a 0 to 1,400-kv 
generator [9]. The potential to the sections was 
graded uniformly by a between 


series of resistors 


the leads from the generator. The 90° beam of 
X-rays (fig. 2) passed through the wall, about 
3-mm copper) of the X-ray tube, by a remotely 
controlled shutter, 5, through a collimator in front 


of the parallel-plate ionization chamber, and through 


& monitor ionization chamber, M. A _ focal-spot 
size of about 2-cm diameter was obtained in the 
X-ray tube. Although this spot size was larger 
than desired, and sometimes shifted during the 
run, the monitor provided a means of normalizing 
the parallel-plate chamber readings because it 


measured most of the beam “seen” by the parallel- 
plate chamber. 

A lead-lined box contained the plates of the ioni- 
zation chamber and limited the X-radiation measured 
by the chamber to that coming through the aperture. 
The box size was determined by pessimistic extr: apo- 
lation from the requirements for lower 


The lead protection required in the box was com- 


energies 





puted according to the method used in NBS Hand. 
book 50 [10] so that the ionization from the lea 
radiation would be than 0.1 percent of the 


less 


ionization produced by the collimated beam. The 
front face of the chamber was shielded by 1... 
thick lead; the sides, top, and bottom by % in. 
thick lead; and the back by \-in.-thick lead. The 


aperture was cut in a 2-in.-thick lead dingheaal 
A 2-in.-thick lead plug over the aperture of the 
diaphragm (with the X-ray beam on) experimentally 
confirmed the adequacy of this protection. The 
snout on the box was provided to reduce the size 
of the box because it was anticipated that a scatter. 
ing diaphragm would be necessary in order to stop 
radiation seattered from the primary diaphragm 
from entering the measuring volume. 
subsequent tests indicate, in agreement with tests 
of Attix and DeLaVergne [11] at lower energies, that 
this diaphragm is not needed. 

The high-voltage plate was a single sheet of alum. 
inum. The collector and field guards were fashioned 
from a single sheet (90 cm by 90 cm) of colloidal 
graphite-coated polyethylene. A collector (30 em 
long by 80 em high) was defined by lines scratched 
through the graphite with a sharp stylus. This 


collector was divided into seven 5-cm strips on each | 


side of a central 10-cm strip by lines scratehed 
through the graphite (see fig. 2). Any or all of 


them could be connected to the ionization-current- 
measuring instrument. When _ not connected, 
they served as part of the guard system. There 
was a 5-cm guard strip above and below the collector 
strip, and a 30-cm-long guard strip in front and in 
back of the collector strips. The plate separation 
could be varied from about 15 to 90 em, and the 
dimensions of the box permitted a movement of 
the plate system over a distance of about 1 m along 
the beam. 

The distorting effect of the grounded box on the 
electric field was reduced by a system of guard strips 
and wires. As pointed out previously [8], the strips 
are preferable but wires are required where photon 
absorption is important. Double wires fastened to 
the inner and outer edges of the st rips were therefore 
installed over an area of 50 by 50 cm at the front and 
back of the chamber system (see fig. 3). Strips 
mm by 3% em) were used over the remaining area 
and the top and bottom. With this system the 
grounded box produced less than 0.2. percent dis- 
tortion * when the plate system was located in the 
center of the box. To reduce the field distortion 
when the plate system was used near the front or 
back of the box, the double wires at both ends of the 
plate system were temporarily replaced by strips 
except for a central area of 15 by 50 em. 

A 10-em-i.d., 1%-mm_ thick, colloidal-graphite- 
coated, Bakelite tube was inserted between the plates 
to eliminate collection of the tonization produced 
by the primary photons. This tube was coaxial with 
the beam and extended in front of and behind the 
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tter- Ficure 2. Elevation view of experimental arrangement. 
Stop [, X-ray target; 8, shutter; D, diaphragm; C, collectors; M, monitor 
‘agm 
ver, | plate system. Scribed lines through the colloidal 
Lests graphite along generatrices of the tube insulated the 
that portions of the tube from each other and permitted 
them to be electrically connected to the guard wires. 
um- With this arrangement, the collected ionization is due 
oned only to scattered photons. The tube undoubtedly 
idal- distorted the electric field in its vicinity but the total 
em scattered photon contribution was small (less than 
rhed | percent of the electron contribution from primary 
This photons), so this distortion was assumed to be negli- 
pach gible. 
‘hed A separate arrangement, similar to that used by 
I of Failla [12] to investigate the ‘electron cloud,” was 
ent- used to measure the growth of electronic equilibrium 
ted, along the beam. A grid chamber of colloidal- 
here graphite-coated nylon thread (~0.2-mm diameter) 
ctor forming 5 parallel, 90 by 90-cm planes was placed 
d in | in the shielded box. The thread spacing was about 
tion | em in each plane, and the planes were spaced about 
the every 1.6 em and perpendicular to the beam. The 
t of center plane (fig. 4) was graphite-coated to form two 
long collectors insulated from each other and from ground 
Piowre 3. Part of free-air chamber showing plate system. | hy uncoated nylon thread. One of these collectors 
the was slightly larger than the photon beam and con- 
rips centric with it. The other included the rest of the 
rips plane. The planes immediately on each side of. this 
ton center plane were connected to the collecting poten- 
| to tial (90 v). The outer two planes of thread were 
lor connected to ground and used to shield the collecting 
and system from the grounded box. 
’? + + One sheet of 1%-mm-thick 90-cm-square conduct- 
_— ing plastic was placed in front of and parallel to the 
the X-RAY plane of the grid chamber. A second was placed 
dis- DIAMETER, | 5 behind the chamber. Each sheet had a center hole 
th MPN slightly larger than the photon beam and each could 
10 s! be moved along the beam. As each of these sheets 
val ” was thicker than the electron range, the ionization 
ips | Seeorageren oy electrons can te ole it thes anaul SaEEE 
SROUND . pas ms : Mee, os 
; " sheets. The growth of ionization was explored by 
ite changing the location of the sheets. 
ites < 
- 3. Procedure 
il 
th = In order to eliminate the effects of X-ray output 
Fictre 4. Schematic elevation of the center plane of the thread | Variation, simultaneous ionization measurements 
chamb were made with the parallel-plate chamber and the 





hit line nd teu ted nylon thread ‘ The large size of the hole makes the exact position somewhat ambiguous. 
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ELECTRON CONTRIBUTION, ARBITRARY UNITS 
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PLATE SEPARATION OR HEIGHT, cm 
Figure 5. lonization collected in air from 400-kv X-rays as 
a function of plate separation with a fixed height, and of 


height with a fixed plate separation. 


, Height; @, plate separation. 


monitor. The charge from each was collected on a 
capacitor and the resulting potential measured under 
a null condition by a vibrating-reed electrometer. 
The ratio of the voltage on the pi arallel- pli ite capac i- 
tor to that on the monitor capacitor was proportional 
to the parallel-plate current for a fixed amount of 
monitor ionization. This averaging of the parallel- 
plate current and the normalization by means of the 
monitor overcame the undesirable fluctuations in the 
X-ray output. Data were obtained at each 50 kv 
in the 250- to 500-kv range. 

The current collected by the collector plate is 
made up of two components, that part due to the 
ionization from high-speed electrons that are pro- 
duced by the primary photons (primary contribu- 
tion), and that part due to ionization from high- 
speed electrons produced by scattered and fluores- 
cent photons (secondary contribution). The sum 
of the two contributions - the secondary contribu- 
tion alone are measured [2 13] so that the primary 
contribution can be de Aa oe by subtraction. In 
the present experiment, this sum was determined 
both by measuring the ionization collected by all of 
the collector strips (80 em high) for different plate 
separations, and by measuring the current collected 
by each strip for a fixed 80-cm_ plate separation. 
The ionization collected by all of the strips at a given 
plate separation, a, should be equal to the ionization 
collected at 80-cm plate separation by all of the 
strips up to a height equal to a. Figure 5 shows a 
representative comparison of the results by the two 
methods. Generally, the agreement was within 0.5 
percent. However, the ionization density decreases 
rapidly with distance from the photon beam, so the 
individual strip method gives a more accurate 
measure of this density function. Therefore, the 
remainder of the data presented here were obtained 








by the strip method. Similar measurements were 
made with the Bakelite tube inserted to obtain the 
secondary contribution. 

The current to each of the collectors of the grid 
chamber was measured separately. The one not 
used for collecting ionization was grounded. The 
variation of the ionization was determined ag 4 
function of the position of the conducting plastic 
sheets with respect to the grid chamber. 

All measurements were obtained at room tempera. 
ture and pressure, which averaged about 26° C and 
760 mm Hg. 


4. Results 


Figure 6 shows the primary contribution for 5-em- 
high layers of air at different vertical distances from 
the beam. The fractional loss in primary photon 
ionization (electron loss) at a particular height js 
obtained by summing the contributions from infinite 
plate height (obtained by linear extrapolation of the 
curves of fig. 6) to the plate height in question and 
dividing by the sum of all of the contributions out 
to infinity. 

The diameter of the photon beam obviously affects 
the loss for a given plate separation. In figure 7 the 
percentage of electron loss is plotted against the 
plate separation minus the penumbra diameter 
(6.3 em) at the center of the collecting region. On 
the 250-kv curve are plotted data taken from results 
of Attix [11],° where the penumbra was about 1.8 


em in diameter and the radiation filtration was 
approximately 1 mm of copper plus 4 mm of alumi- 
num. This approximate method of treating the 


beam diameter seems to be adequate at least for 
these data, as the agreement of the electron contri- 
butions for the two experiments is generally within 
0.1 percent. Furthermore, the difference is in the 
expected direction because the higher filtration used 
in the present experiment results in a smaller number 
of high-energy photoelectrons. 

In order to compare the present data with those 
published earlier by Block [13] and to provide more 
vasily usable data for other chamber configurations 
than the parallel-plate type, the electron loss and 
the scattered photon contribution were computed 
for cylindrical coordinates. Obviously the ioniza- 
tion collected from a ring of a given radius will be 
independent of angle, so only one variable is re- 
quired. The contribution to each 5-em ring was 
obtained from the results shown in figure 6 by the 
method previously outlined [4]. Figure 8 shows the 
electron loss per 5-cm-thick ring obtained at different 
inner radii. Here again the losses may be summed 
from infinity to a given radius to obtain the loss 
beyond that radius. The results are shown iD 
figure 9, where the beam size has been reduced to 
zero by subtracting the beam radius. 

Figure 9 also shows some points computed from 
the data of Block [13] at 400 kv and a filtration of 


These pl d data are slightly different from those in table 2 of reference [4 


because the nlotted data ire for an effectively infinite plate heig 
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Figure 8. Percentage primary electron tonization loss 


per 5-cm ring for different inner radii from the beam. 


of 


The beam penumbra diameter was 6.3 em at the center of the collecting region. 
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Results of Attix for 250 kv 


2.55 mm of tin plus 0.5 mm of copper plus 5.0 mm 
of aluminum. These points were obtained with a 
different filtration than that used for the present 
data; some of these had to be determined from small- 
scale drawings, the scattered photon contribution 
had to be computed from the present data, and the 
beam diameters were assumed to be the same for 
Block as for the present work. It is thus not sur- 
prising that these points differ by up to 0.5 percent 
from the present work. 
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Figure 9. Percentage of primary 


beyond different radii from the 


Results corrected to zero beam diameter , Data of Attix [4] at 250 kv; A, 


data of Block [13] at 400 kv 


Figure 10 gives the radial distribution of the ioniza- 
tion from the secondary photons obtained from 
measurements made with the Bakelite tube in a 
fashion similar to that for the primary contribution. 
The contribution is given as a percentage of the total 
ionization from primary electrons. As the curves 
for various X-ray-tube potentials have a_ similar 
shape, only the 250-kv curve is shown. The con- 
tribution at 300-kv is obtained by multiplying the 
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\ -ray be an 
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, Present « 


I 100 kv, multi; the ordinate by 
0.9: for 400 and 


ita , data of Attix [11 For 
S00 kv, multiply by 0.8 


ordinates of the curve by 0.9; at 400 and 500 kv, by 
0.8. It is interesting to note that the data obtained 
by Attix [11] at 250-kv on a much smaller chamber 
give very nearly the same results. Here again the 
abscissa has been corrected for beam size by sub- 
tracting the beam radius. 

The following table shows the air absorption per 
meter at various tube potentials obtained from ioniza- 


tion-chamber measurements with the chamber at 
the front and at the back of the box. Repeat 
measurements agreed to within 0.2 percent. 


They also agreed to within 0.2 percent with the results 


of Block 


Tube Nir 


potential absorption 


Ai ” 
250 | 
350 > 
500 13 


Figure 11 shows a typical “growth-of-electronic- 
equilibrium” curve for the ionization collected by 
the outer electrode of the thread chamber with 
500-kv X-rays. The charge collected by the center 
electrode is 0.153 arbitrary unit, and is independent 
of the position of the plastic sheets. When the 
sheets are at a large distance from the collecting 
electrodes the total (equilibrium) ionization collected 
is 0.0444+-0.032—0.153— 0.0296, or about 0.200 
arbitrary unit. From these data for 500 kv, it is 
seen that (0.153 0.029)/0.200. or about 90 percent, 
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of this ionization results from electrons effectively 
originating within a few centimeters of the collecting 
region. Approximately 1 percent of the equilibrium 
ionization is produced by electrons originating up to 
about 15 em behind, and about 90 percent by elec- 
trons originating up to 40 em in front of the collecting 
region. The flux of gamma rays producing these 
high-speed electrons of course varies over the dis- 
tance (404-15 em) because of air absorption. 

To a first approximation the absorption of X-rays 
in the air between the diaphragm and the origin of 
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the high-speed electrons may be computed from the 
gir-absorption coefficient and the distance between 
the diaphragm and the collecting region. A correc- 
tion to this distance may be obtained from curves 
such as those in figure 11. For 500-kv X-rays the 


correction is 
1 "40 

| +-(0.0444)(40—3.7)— | li | 
0.200 J3.7 


715 
+ (0.032 (15 3.7)— | rai | 


where 0.200 is the equilibrium ionization and di is 
an increment of ionization at a distance z, in front of 
the collector and x, behind the collector. On substi- 
tuting values from the curve, this distance correc- 
tion is found to be less than 1 em. As 8 em of air 
attenuates the photon beam by only the order of 
0.1 percent, the actual displacement is negligible. 
For all lower X-ray potentials the displacement is 
smaller. Therefore, the air-absorption correction 
in the 250- to 500-kv range should usually be com- 
puted from the center of the collector to the aperture. 

Figure 12 shows the approximate distance required 
for electronic equilibrium. Also plotted are the dis- 
tances equal to the radii from the beam, at which the 
loss of ionization from the primary electrons is just 
compensated by the gain of ionization from the 
secondary radiation. Thus, a collector-to-aperture 
distance, equal to the radius for compensation, is ade- 
quate for moderately filtered X-rays of 400 kv and 
below but may be insufficient for 450 and 500 kv. 
It is seen from figure 11, for instance, that a decrease 
n the collector-to-aperture distance from 40 to 29.5 
em results in an error of about (0.0444—0.0438 
0.200, or 0.5 percent. 
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5. Summary 


The data required for the design of free-air stand- 
ard chambers for measurement of 250- to 500-kv X- 
rays in roentgens have been obtained. The internal 
consistency of these data indicates that the errors of 
plate spacing can be estimated to about 0.2 percent. 
However, the much closer agreement between these 
results and those of Attix indicates that the uncer- 
tainty is probably less than 0.2 percent. 
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